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Section  I 


INTRODUCTION 


Metal-forming  is  a  process  in  which  the  shape  of  a  workpiece, 
consisting  of  a  metal  or  its  alloy,  is  altered  into  useful  shapes,  by 
using  shaped  dies  under  high  pressure.  The  knowledge  of  the  detailed 
deformation  mechanics  under  various  working  conditions  and  the  accurate 
determination  of  the  forming  load,  as  well  as  the  values  of  stress  and 
strain  distributions,  are  extremely  important  to  the  proper  design  and 
control  of  forming  processes.  Because  of  large  amounts  of  plastic  defor¬ 
mation,  the  complexities  in  geometric  changes,  and  the  effects  of  the 
frictional  force  along  the  tool -workpiece  interfaces,  the  exact  mathe¬ 
matical  closed-form  solution  for  the  metal-forming  problem  is  usually 
very  difficult  to  obtain,  except  for  some  special  cases.  Most  of  the 
closed-form  solutions  are  therefore  based  on  many  assumptions,  which 
may  oversimplify  the  actual  situations  and  limit  their  value  in  practical 
applications . 

The  closed- form  solutions  were  usually  made  by  one  of  the  following 
methods  [1],  [2],*  namely,  Sach's  [3],  or  slab,  method;  uniform  deforma¬ 
tion  energy  method  [1],  [4];  slip-line  method  fl],  [2],  [5],  [6];  and 
bounding  method  [1],  [2],  [7],  [8].  In  the  slab  method,  the  homogeneous 
deformation,  which  states  that  the  stresses  only  change  in  one  direction 
and  that  the  frictional  constraint  at  the  tool -workpiece  interface  does 
not  affect  the  internal  stress  distribution,  is  assumed.  The  uniform 


‘Numbers  in  brackets  correspond  to  the  references  listed  at  the  end  of 
this  report. 
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deformation  energy  method  calculates  the  average  forming  stress  from  plastic 
deformation  by  neglecting  the  redundant  work  involved  in  internal  shearing 
processes  due  to  nonuniform  deformations.  Both  of  these  two  methods  can 
only  provide  the  overall  quantity,  such  as  the  forming  load,  but  not  the 
detailed  stress  and  strain  distributions.  Furthermore,  the  assumption  of 
homogeneous  deformation  oversimplifies  the  actual  deformation  pattern  in 
most  forming  processes. 

The  first  approach  to  the  analysis  of  metal-working  processes  that 
did  not  assume  homogeneous  deformation  was  the  slip-line  method.  This 
method  permits  a  point-by-point  calculation  of  stress  for  problems  of 
plane-strain  deformation  in  an  idealized  non-work-hardening  rigid-plastic 
solid.  Although  the  method  provides  information  about  local  stress  and 
velocity  distributions,  the  method  itself  lacks  uniqueness  for  constructing 
slip-line  fields.  Furthermore,  the  solutions  are  limited  to  rigid- 
perfectly  plastic  material  under  plane-strain  forming  conditions. 

The  bounding  methods  give  estimated  upper  and  lower  limits  of  the 
required  deformation  force,  according  to  kinematically  admissible  velocity 
fields  (upper-bound)  or  statically  admissible  stress  fields  (lower-bound) 
provided  for  the  forming  processes.  In  order  to  have  a  mathematical  closed- 
form  solution,  the  provided  admissible  field  must  be  mathematically  expres¬ 
sible  and  the  rigid-perfectly  plastic  material  behavior  is  usually  assumed. 
The  construction  of  these  admissible  fields  is  usually  very  difficult, 
especially  for  geometrically  complex  forming  processes.  In  addition,  these 
methods  lack  the  ability  to  efficiently  deal  with  work-hardening  materials 
and  to  reveal  detailed  deformation  mechanics,  such  as  the  local  stress  and 
strain  distributions. 
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In  the  last  several  years,  due  to  rapid  developments  in  computers, 
the  closed-form  solutions  have  become  unnecessary,  and  various  numerical 
techniques  have  been  used  for  the  analysis  of  metal-forming  processes. 

The  well-known  finite-element  methods  [9],  [10],  [11]  were  applied  to 
the  problems  of  deformation  in  elastic-plastic  solids  for  such  processes 
as  flat-punch  indentation,  upsetting,  side  pressing,  and  extrusion  [12] - 
[16].  The  elastic-plastic  finite-element  method  offers  many  advantages  in 
providing  detailed  deformation  mechanics  and  in  obtaining  accurate  solutions. 
Work-hardening  material  behavior  and  the  interface  friction  constraint  can 
also  be  treated  efficiently  by  this  method.  In  spite  of  the  many  advantages 
provided  by  this  method,  however,  there  s'till  remain  serious  restrictions 
regarding  computations.  Since  the  method  deals  with  elastic-plastic 
materials,  elastic  calculations  are  ne  ssary  at  every  step,  as  is  constant 
checking  to  determine  if  any  element,  previously  elastic,  has  reached  the 
yielding  stage.  When  an  element  reaches  the  yielding  stage,  the  stiffness 
relations  for  that  element  change,  thereby  limiting  increment  size. 
Furthermore,  the  elastic-plastic  property  itself  does  not  permit  the  use 
of  large  increments  of  deformation.  Because  of  this  requirement  of  small 
increments,  the  extensive  computation  time  necessary  to  apply  this  method 
to  metal-forming  problems  renders  it  cumbersome  and  uneconomical. 

Recently,  based  on  a  variational  principle  similar  to  the  upper-bound 
approach  for  rigid-plastic  materials,  Lee  and  Kobayashi  [17]  have  devel¬ 
oped  a  method  called  the  rigid-plastic  finite-element  method  (also  called 
the  matrix  method).  This  method  not  only  has  all  the  advantages  of  the 
finite-element  method,  but  it  also  permits  large  increments  during  deforma¬ 
tion,  therefore  reducing  computation  time.  Successful  applications  of 
this  method  were  made,  for  example,  in  the  prediction  of  workability  in 
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upsetting  [18],  [19]  and  the  analysis  of  plane  stress  bore  expanding  [20]. 
Modifications  of  this  method  were  made  for  sheet-metal  problems  in  which 
large  deformation  and  rotation  are  involved  and  the  geometric  nonlinearity 
cannot  be  neglected  [21].  Parallel  to  the  rigid-plastic  finite-element 
method,  based  on  a  similar  variational  principle,  Godbole  and  Zienkiewicz 
[25]- [27]  and  Price  and  Alexander  [22]- [24]  were  able  to  solve  the  rigid- 
viscoplastic  and  viscous  flow  problems  of  plane-strain  extrusion,  compres¬ 
sion,  and  flat -punch  indentation. 

To  successfully  apply  the  rigid-plastic  finite-element  method  to  the 
analyses  of  various  forming  processes,  it  requires  suitable  modifications 
according  to  the  individual  problem  characteristics.  An  efficient  method 
for  dealing  with  the  rigid  elements  inside  the  control  volume  is  also 
required  in  order  to  avoid  bad  convergence  characteristics.  With  a 
suitable  choice  of  the  convergence  criterion,  however,  the  computing  time 
required  to  obtain  the  convergence  solutions  can  also  be  reduced  signifi¬ 
cantly.  The  present  investigation  is  therefore  aimed  at  improving  and 
modifying  this  method  to  analyze  various  types  of  problems,  such  as 
steady-state  problems  of  extrusion  and  drawing,  neutral  flow  problems  of 
ring  compression,  and  the  problems  of  complex  geometric  change  in  closed- 
die  forging. 

In  this  report.  Section  II  discusses  the  minimum,  or  stationary 
value,  principle  and  the  methods  of  imbedding  the  incompressibility 
constraints  for  rigid-plastic  materials.  The  boundary  conditions  along 
the  die-workpiece  interface  and  the  methods  of  rigid-body  treatment  are 
also  discussed.  Section  III  describes  the  finite-element  method  and 
Newton-Raphson's  iterative  technique  used  in  obtaining  the  minimum,  or 
stationary  value,  solutions.  Section  IV  investigates  the  deformation 
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mechanics  of  the  multipass  steady-state  extrusion  and  drawing  processes, 
emphasizing  the  influence  of  the  previous  deformation  history  on  the  cur¬ 
rent  metal  flow. 

In  implementing  the  finite-element  method  for  the  analysis  of  metal¬ 
working  problems,  particular  attention  must  be  paid  to  boundary  conditions. 
In  upsetting,  a  part  of  an  initially  free  surface  comes  into  contact  with 
the  die  during  deformation.  A  similar  situation  occurs  in  sheet-metal 
forming.  There  is  also  a  class  of  problems,  such  as  ring  compression  and 
rolling.  A  unique  feature  of  these  kinds  of  problems  is  that  there  exists 
a  point  (or  a  region)  along  the  die-workpiece  interface  where  the  velocity 
of  the  deforming  material  relative  to  the  die  velocity  becomes  zero  (the 
problems  of  neutral  flow  point  or  region).  The  frictional  stress  usually 
changes  its  direction  at  the  neutral  flow  point,  but  the  lcoation  of  this 
point  is  not  known  a  priori.  In  Section  V,  a  variational  formulation  is 
developed  for  the  case  where  the  frictional  stress  is  dependent  on  the 
relative  velocity  at  the  die-workpiece  interface.  This  modified  formula¬ 
tion  is  then  used  to  analyze  the  ring  compression  problems. 

The  closed-die  forging  process  is  an  extremely  difficult  problem 
because  of  the  complex  die  geometries,  the  large  rigid  portions,  and  the 
rapid  changes  in  velocity  field  during  deformations.  In  Section  VI,  the 
rigid-plastic  finite-element  method  is  applied,  with  further  modifications, 
to  analyze  these  extremely  complex  problems. 
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Section  II 


PROBLEM  FORMULATION 


I.  MINIMUM  PRINCIPLE 

The  Cartesian  coordinate  is  used  to  derive  the  minimum  principle. 
Generalization  of  this  formulation  to  other  coordinate  systems  can  be 
made  without  difficulty.  A  body  of  volume,  V,  is  considered  with  the 
traction,  F,  prescribed  on  a  part  of  the  surface,  Sp,  and  the  velocity, 

V,  prescribed  on  the  remainder  of  the  surface,  S^.  The  body  is  composed 
of  a  rigid-plastic  material  which  obeys  the  von  Mises  yield  criterion 
and  its  associated  flow  rule.  Body  forces  are  assumed  to  be  absent.  The 
actual  stress  (a^)  and  velocity  field  (v)  satisfy  the  following  relations 

(a)  Equilibrium  conditions: 


0.  .  .  =  0; 
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(b)  Compatibility  conditions: 


e. .  =  T(v.  .  +  v.  . ) 
lj  2  i,j  j,i 


and 


(1) 


(2) 


e  =  v.  .  =  0; 
v  1,1  ’ 


(c)  Stress-strain  rate  relationship: 


o! .  =  4  t  e.  •; 
13  3-  13 


(3) 
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where  o! .  are  the  deviatoric  stresses,  0  and  I  are  the  effective 
ij 

stress  and  the  effective  strain  rates,  respectively,  and  are 

defined  by  o  =  J ^  a ! . a ! .  and  e  =  J \  c.  .1-  ■ . 

2  1J  ij  3  i] 

(d)  Boundary  conditions: 


a.  .n. 
ij  i 


v. 

l 


on 


on 


S 

S 


F 

U 


where  p  is  the  unit  outward  normal  to  the  corresponding  surface. 


Let  v*  be  any  admissible  velocity  field .satisfying  the  volumetric  con- 
straint  e*  =  0  and  the  prescribed  velocity  boundary  conditions.  Multiplica¬ 
tion  of  the  left-hand  side  of  Eq.  (1)  by  v?  and  integration  of  this  quantity 
over  the  whole  volume  gives  the  following  virtual  work  principle: 


L  °iiha  dv '  L  Vi ds 


*  . 


(a.  . n . ) V.  dS, 
13  J  i 


(4) 


where  eV  are  the  strain  rate  components  compatible  with  the  admissible 
velocity  field  v*.  Note  that  in  deriving  Eq.  (41  the  divergence  theorem 
is  used. 

The  convexity  and  normality  conditions  of  the  yield  surface  provide 
the  following  inequality. 


(cij  -  i  o,  (1 

where  are  the  stress  components  derived  from  the  admissible  velocity 
field  v*  through  the  stress-strain  rate  relationship  (3),  satisfying  the 
yield  criterion. 
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Because  the  volumetric  constraint  e*  =  0  is  imposed  on  the  admissible 
velocity  field,  the  tensor  product  of  the  stress  and  strain  rate  is  there¬ 
fore  equivalent  to  the  product  of  effective  stress  and  effective  strain 
rate  and  is  expressed  by 


0*.  .  £  *  .  -  ci* 
11  1J 


(6) 


Substitution  of  Eqs.  (6)  and  (5)  into  Eq.  (4)  gives  the  following 
inequality: 

f  5  e*  dV  -  f  F*v*  dS  > 

’F 


ff  £*  dV  -  J  F*v*  dS 
V  S. 


(c. .n.JV.  dS, 
nil 


(7) 


where  the  equality  holds  only  if  ^ ,  i.e.,  when  the  stress  compo¬ 

nents  correspond  to  the  actual  solutions.  The  functional  4>  defined  by 


*  ■  L 


c  e*  dV  - 


J, 


F  •  v  *  dS 


(8) 


is  therefore  used  for  the  finite-element  formulation,  which  states  that, 
among  all  admissible  velocity  fields  v*,  the  actual  solution  gives  the 
minimum  value  of  this  functional  <J). 

This  minimum  principle  is  also  found  in  the  books  by  Hill  [5]  and 
Washizu  [28]  and  has  been  used  extensively  in  the  analysis  of  metal -forming 
processes  [17],  [29].  As  pointed  out  by  Washizu  [28]  and  Kobayashi,  Lee, 
and  Oh  [30],  the  solution  obtained  from  the  minimization  of  the 

functional  (8)  differs  from  the  actual  strain  rate  by  a  scalar  factor 

•  • 

c  (i.e.,  =  ce^j),  if  surface  tractions  are  prescribed  over  the  entire 

surface:  and  the  solution  a*,  differs  from  the  actual  stress  o. .  by  a 

ij  lj 

uniform  hydrostatic  stress  p  (i.e.,  +  p6„),  if  velocities  are 
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prescribed  over  the  entire  surface,  or  if  a  combination  of  tangential 
traction  and  normal  velocity  is  prescribed  over  that  surface. 

Although  the  formulation  above  is  still  an  upper-bound  approach,  the 
method  used  in  obtaining  the  solution  is  quite  different.  The  finite- 
element  techniques  used  in  the  present  analysis  gives  much  more  attractive 
results  than  those  of  classical  upper-bound  solutions.  In  the  classical 
upper-bound  predictions,  an  admissible  velocity  field  is  provided  first. 
This  admissible  field  must  be  mathematically  expressible  in  order  to  have 
a  mathematical  closed-form  solution.  Although  the  provided  velocity  field 
may  be  verified  experimentally,  a  reasonable  mathematical  representation 
of  this  velocity  field  is  usually  very  difficult  to  obtain.  In  addition, 
the  classical  upper-bound  solution  lacks  the  ability  in  efficiently  deal¬ 
ing  with  work-hardening  materials;  it  also  lacks  the  abilities  in  reveal¬ 
ing  the  detailed  deformation  mechanics,  such  as  the  local  stress  and  strain, 
and  the  die  pressure  distribution.  These  severe  restrictions  limit  the 
value  in  the  actual  application  of  the  classical  upper-bound  method.  On 
the  contrary,  the  present  analysis  makes  use  of  the  finite-element  tech¬ 
niques  to  solve  the  functional  (8)  above;  the  volumetric  constraint  is 
removed  by  introducing  the  Lagrange  multiplier  or  penalty  function  formula¬ 
tions  (see  Section  I I. 4).  Certain  type  of  velocity  distributions  inside 
an  element,  in  terms  of  surrounding  nodal  velocities,  is  presumed;  but 
the  nodal  velocities  themselves  are  obtained  from  the  computational  results 
directly.  Work -hardening  material  behavior  and  frictional  force  between 
the  die-workpiece  interface  can  also  be  treated  efficiently.  Besides, 
this  method  also  provices  the  details  necessary  for  the  proper  design  and 
control  of  the  metal-working  process. 


2.  RIG ID- BODY  TREATMENT 


The  minimum  principle  described  in  the  previous  section  applies  only 
if  the  entire  body  is  deforming  plastically  and  no  rigid  zone  or  unloading 
exists  during  the  deformation  process.  In  practical  problems,  however, 
situations  do  arise  where  the  rigid  zone,  as  well  as  rigid  unloading,  are 
involved.  These  rigid  portions  are  characterized  by  the  extremely  small 
value  of  effective  strain  rate  in  comparison  to  that  in  the  deforming  body. 
If  these  portions  are  included  within  the  control  volume,  V,  the  value 
of  the  first  derivatives  of  the  functional  (8)  cannot  be  uniquely  deter¬ 
mined  because  of  the  undefined  value  of  the  effective  stress  (Eq.  (3)) 
when  the  effective  strain  rate  approaches  zero.  Due  to  these  undetermined 
derivatives,  solutions  are  impossible  to  obtain,  except  for  the  case  where 
the  rigid  portions  are  eliminated  from  the  actual  calculation.  However, 
the  rigid  elimination  technique  is  virtually  inapplicable  in  most  of  the 
problems  because  of  the  difficulties  in  determining  the  rigid-plastic 
boundaries . 

In  an  effort  to  eliminate  the  problem  of  the  undetermined  functional 
derivatives,  the  following  method  makes  use  of  a  limited  value  (say, 
eQ  =  10  )  in  Eq.  (3)  for  all  elements  lying  inside  the  rigid  portions. 

Elements  with  the  value  of  effective  strain  rate,  e,  less  than  the  as- 
signed  value,  e^,  are  called  near-rigid  elements.  These  elements  are 
not  eliminated  from  the  control  volume  during  actual  computations.  The 
limited  value,  eQ,  assumed  in  the  stress-strain  relationship  for  these 
elements,  is  to  ensure  a  unique  value  of  the  deviatoric  stress.  This 
presumptive  stress-strain  rate  relationship  is  equivalent  to  the  assumption 
of  a  Newtonian  fluid-like  material  behavior  for  those  elements  lying  inside 
the  near-rigid  portions  and  is  expressed  by 
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first  derivatives  with  respect  to  v*,  and  the  second  derivatives  are  only 
constants  depending  on  a  and  E^. 

This  modified  functional  (12)  was  used  in  the  present  study.  In 

problems  such  as  extrusion  and  drawing,  the  near-rigid  portion  does  not 

present  much  difficulty  in  obtaining  the  convergent  solutions,  since  the 

boundary  conditions  are  mainly  velocities  prescribed.  For  the  forging 

problems,  however,  a  large  portion  of  force-free  boundary  is  prescribed. 

This  force-prescribed  boundary,  together  with  the  large  rigid  portions, 

results  in  bad  converging  characteristics.  The  present  method  resolves 

this  difficulty  by  means  of  a  suitable  selection  of  the  E^  value  (say, 

-  -4  -  3 

Gp  =  10  in  extrusion  and  10  in  forging)  and  results  in  reliable  solu¬ 
tions  with  fast  convergence. 

3.  BOUNDARY  CONDITIONS  ALONG  THE  TOOL-WORKPIECE  INTERFACE 

In  general,  neither  velocity  nor  force  (including  magnitude  and  dir¬ 
ection)  is  prescribed  along  the  tool -workpiece  interface.  This  is  because 
the  direction  of  the  frictional  stress  is  opposite  to  the  direction  of 
the  relative  velocity  between  the  workpiece  and  the  die,  and  this  relative 
velocity  is  not  known  a  priori.  However,  situations  exist,  e.g.,  in  ex¬ 
trusion  and  drawing,  in  which  the  direction  of  metal  flow  is  well  under¬ 
stood.  This  gives,  in  turn,  the  direction  of  the  frictional  stress.  This 
class  of  problems  can  be  solved  if  the  magnitude  of  the  frictional  stress 
is  given  according  to  the  well-known  Columbus  Law,  expressed  by 

f  =  MP,  ‘  (12) 

or  the  friction  law  of  constant  factor  m,  expressed  by 

12 
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f  =  mk. 


(13) 


where  p  is  the  die  pressure  and  k  is  the  shear  yield  stress.  The  tool- 
workpiece  interface  is  treated  as  a  force-prescribed  boundary  and  the 
formulations  in  the  sections  above  are  readily  applicable.  For  metal¬ 
forming  problems,  the  friction  law  of  Eq.  (13)  has  been  used  extensively; 
examples  are  found  in  [29],  [39],  and  [47], 

For  problems  such  as  ring  compression  and  rolling,  the  unknown  direc¬ 
tion  of  the  relative  velocity  between  the  die-workpiece  interface  makes 
it  extremely  difficult  to  handle  the  boundary  condition.  A  unique  feature 
of  this  type  of  problem  is  that  there  exists  a  point  (or  a  region)  along 
the  die-workpiece  interface  where  the  velocity  of  the  deforming  material 
relative  to  the  die  becomes  zero  (a  neutral  flow  problem),  and  the  loca¬ 
tion  of  this  point  (or  region)  depends  on  the  magnitude  of  the  frictional 
stress  itself.  To  deal  with  these  problems,  a  new  variational  formula¬ 
tion  has  been  developed  and  successfully  applied  to  the  case  of  ring  com¬ 
pression.  The  details  of  these  new  derivations  will  be  discussed  in 
Section  V. 

4.  IMBEDDING  OF  THE  INCOMPRESSIBILITY  CONSTRAINT 

The  incompressibility  condition,  e*  =  0,  imposed  on  the  velocity 
fields  can  be  removed  by  introducing  either  the  Lagrange  multiplier  [17], 
[30],  or  a  penalty  function  [31]  into  the  functional  <t>  (Eq.  (8)).  Note 
that  the  minimum  principle  becomes  a  stationary  value  problem  if  the 
Lagrange  multiplier  method  is  used.  In  order  to  avoid  the  confusion  of 
the  notations,  the  variable  u  is  assigned  to  represent  any  set  of  the 
arbitrary  velocity  fields  which  satisfies  the  velocity  boundary  conditions 
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Comparing  the  variations  of  the  functionals  and  the  penalty  func¬ 
tion  approach  can  also  provide  the  mean  stress,  if  the  velocity  u  obtained 
from  the  minimum  value  of  is  identical  to  that  obtained  from  the 
stationary  value  of  and  this  mean  stress  is  given  by 

dv  •  A  =  o  .  (17) 

v  m 

Both  approaches  have  advantages  and  disadvantages,  as  follows:  (i)  For 
the  Lagrange  method,  the  additional  unknown  A  is  introduced  and,  if  a 
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large  number  of  elements  is  used  in  the  finite-element  discretization, 
this  additional  unknown  will  result  in  a  large  bandwidth  for  the  stiff¬ 
ness  matrix  and  also  requires  considerable  amounts  of  computing  time  to 
solve  the  equation;  (ii)  For  the  penalty  approach,  although  the  additional 
unknown  is  avoided  in  the  problem  formulation,  the  system  of  equations 
may  become  ill-conditioned  when  a  bad  initial  guess  is  provided,  giving 
a  very  large  value  of  ££v.  It  is  suggested,  therefore,  to  use  the  upper- 
bound  velocity  distribution  as  an  initial  guess,  if  available. 

For  a  four-nodal  quadrilateral  element  which  is  used  in  the  finite- 
element  discretization,  it  is  shown  (see  Section  III. 5)  that  the  volumetric 
strain  rate  is  linearly  distributed  inside  that  element.  It  can  be 
shown  further,  for  the  case  of  plane  stress  or  plane  strain,  that  the  in¬ 
compressibility  constraint,  £  dV  =  0,  can  always  be  achieved  if  the 

Jy  v 

volumetric  strain  rate  is  varying  linearly  from  positive  to  negative  with 
a  zero  value  at  the  center  of  the  element,  without  any  restriction  on  its 

f  *2 

slope.  The  penalty  function  constraint,  e  dV,  on  the  other  hand,  re- 

■'V  v 

quires  the  value  of  £y  to  be  zero  at  every  internal  point  of  the  element 

•2 

because  of  the  square  quantity,  £y,  of  the  integrand.  Consequently,  due 
to  this  over-constraint  requirement,  an  incorrect  solution  may  result,  if 
the  penalty  function  (15)  is  applied  without  suitable  modifications.  The 
well-known  "reduced  integration"  has  been  proved  to  be  successful  in  deal¬ 
ing  with  the  penalty  term  for  the  cases  of  plane  stress  and  plane  strain 
and  for  three-dimensional  problems  [9] .  This  reduced  integration  releases 

the  over-constraint  conditions  by  the  method  of  reducing  the  integration 

•2 

points  in  numerical  evaluation  of  the  integral,  e  dV.  For  example, 

Jy  v 

the  2x2  numerical  integration  is  reduced  to  one-point  numerical  integration 
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for  the  cases  of  two-dimensional  plane  problems;  the  only  condition  re¬ 
quired  to  satisfy  the  volumetric  constraint  is  the  requirement  of  the  zero 

• 

at  the  integration  point  (i.e.,  the  centroid  of  the  element),  but  not 
necessary  at  other  points  (see  Section  III. 5).  It  can  also  be  shown  that 
the  reduced  one-point  integration,  evaluated  at  the  centroid  of  the  element, 


gives  the  exact  evaluation  of  the  integral. 


e  dV  for  cases  of  plane 


stress  and  plane  strain,  but  not  for  axisymmetric  cases.  Nevertheless, 
the  over-constraint  conditions  in  axisymmetric  problems  can  also  be  removed 
by  the  method  of  reduced  integration,  although  the  exact  evaluation  of 


the  integral. 


dV  is  not  attainable  from  this  reduced  integration. 


A  slight  modification  of  the  penalty  term  can  also  provide  the  de¬ 
sired  solutions  without  employing  the  reduced  integration.  This  modified 
penalty  form  is  stated  as 


4>,  =  |  oe  dV  +  ~  [  I  £  dV 

1  LJ  v 


J. 


F*y  dS 


(18) 


and,  when  the  admissible  velocity  field  corresponds  to  the  actual  one, 
the  mean  stress  is  given  by 

‘  ■  1 1  K  «•  (19) 


Note  that,  in  Eq.  (19),  the  mean  stress  is  assumed  to  be  uniform  over 
each  element  as  that  assumed  for  the  Lagrange  multiplier  approach.  Com¬ 
parison  of  Eqs.  (19)  and  (17)  gives  the  facts  that  the  present  modifica¬ 
tion  is  made  by  assuming  the  average  volumetric  strain  rate 


;ave 

C. 

V 


dV 


(20) 
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and  by  substitution  of  this  average  value  into  Eq.  (15).  The  square  out 
side  the  integral  e  dV  makes  the  formulation  free  from  the  possibil- 

Jv  v 

ities  of  over-constraint  and  renders  this  formulation  successful  in  appl 
cation  to  any  problem. 

In  the  present  investigations,  the  Lagrange  multiplier  method  (14) 
is  applied  to  the  problems  of  ring  compression;  the  modified  penalty 
approach  (18)  is  used  in  the  analysis  of  multipass  extrusion  and  drawing 
and  forging  problems.  Further  details  in  terms  of  finite-element  dis¬ 
cretization  are  given  in  Section  III. 


Section  III 


METHOD  OF  SOLUTION- -THE  FINITE-ELEMENT  APPROACH 


1.  METHOD  OF  APPROACH 

General  theories  of  the  finite-element  method  have  been  well  described 
in  many  books  (for  example,  Zienkiewicz  [9],  Strang  and  Fix  [10],  Bathe 
and  Wilson  [32]).  These  finite-element  techniques  were  applied  to  the 
present  minimum  or  stationary  value  problems  for  rigid-plastic  work¬ 
hardening  materials  with  the  functionals  as  given  in  Section  II.  The 
method  is  briefly  described  as  follows. 

Let  the  functional  4>  (w)  represent  a  given  volume  integral  with  its 
value  evaluated  over  the  whole  body,  V.  The  body,  V,  is  divided  into 
M  elements  interconnected  at  N  nodal  points.  The  value  of  the  functional, 
<j>(w),  may  be  written  as 


M 

-  I 

m=l 


*(n)(w). 


(21) 


where  <j>^m^(w)  represents  the  same  integral  but  evaluated  at  the  m-th 

element.  Inside  this  element  the  variables,  w,  are  approximated  by  a 

'  (*) 

linear  combination  of  certain  shape  functions,  N^  such  that 


w  —  J  N'w',  (22) 

L  a  a.' 
a 

where  the  unknown  coefficients  w^  may  be  explained  as  the  values  of  the 

variables  w  at  nodal  point  a,  and  are  called  nodal  variables.  (These 

(*)  Hereafter,  the  Greek  indices  (a,3, . . . ,etc.)  are  reserved  to  denote 
the  surrounding  nodal  points  or  their  associated  variables  at  the 
elemental  level. 
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nodal  unknowns  may  represent  displacements,  velocities,  temperature,  etc., 
depending  on  the  type  of  the  problem).  Substitution  of  Eq.  (22)  into 
Eq.  (21)  gives  the  approximation  of  the  functional  <p  by 


which  depends  on  the  nodal  point  variables  w^  only. 

If  w.  (i  =  1 , . . . , N  x  p;  N  =  total  nodals,  p  =  number  of  unknowns  to 
each  nodal)  denote  the  unknowns  of  the  nodal  variables  (y^)  in  the  global 
system  after  assembly,  the  functional  <p  in  Eq.  (23)  can  then  be  replaced  by 


4>  -  <Kwi)  (i  =  1, . .  .  ,N  x  p) . 


(24) 


Thus,  the  discretization  procedures  change  the  original  stationary  value 
problem  of  the  functional  <p ,  with  respect  to  w  to  the  problem  of  the 
functional  $  with  respect  to  nodal  variables  w^. 

To  find  the  minimum  or  stationary  value  of  the  functional  $  of 
Eq.  (24),  the  following  system  of  equations  must  be  solved: 

94>(w. ) 

— Z~~  =  0  (j  =  1.....N  x  p).  (25) 

9w. 

J 

If  this  system  of  equations  is  linear,  it  is  then  possible  to  write 
Eq.  (25)  in  terms  of  the  following  matrix  representation: 


[K..]{w.}  =  {?.},  (26) 

where  ] ,  known  as  the  stiffness  matrix,  is  independent  of  w^ ;  and 
{f^}  is  the  nodal  force  vector.  The  matrix  ]  is  usually  a  banded 
symmetric  matrix  in  most  solid  mechanics  problems.  For  these  cases, 


only  the  upper  half  of  the  elements  of  matrix  [K„]  are  actually  required 

during  computation,  due  to  the  symmetric  property;  and  these  elements 

can  be  arranged  and  stored  in  a  square  array  inside  the  computer,  due  to 

the  banded  property.  The  Gaussian  elimination  and  backward  substitution 

techniques  are  normally  used  in  solving  the  equations  above  (26) . 

If  the  elements  of  the  stiffness  matrix  [K..]  are  functions  of  the 

i  J 

nodal  variable  w^ ,  the  system  (25)  is  then  called  nonlinear.  All  the 
functionals  given  in  Section  II  belong  to  this  type  of  problem.  Normally, 
it  requires  more  effort  to  solve  the  nonlinear  problems  than  the  linear 
ones.  Most  of  the  methods  rely  on  the  iterative  techniques  [33],  [34]. 
although  the  convergence  character,  during  iterations,  depends  very  much 
on  the  functional  behavior  itself,  the  well-known  Newton-Raphson  iterative 
method  has  been  extensively  used  due  to  its  simplicity  and  the  quadratic 
convergence  character. 

With  the  assumptions  of  the  existence  and  continuity  of  the  derivatives, 
3$ 

the  value  of  -  evaluated  at  an  arbitrary  point  w,  (n)  can  be  written  as 

3w. 

l 

the  following  series  representation  by  virtue  of  Taylor's  expansion: 


3$  !  3$  !  3 “•*>  I 

— —  =  - 1  +  - 1  £w4  ,  . 

3w.'w,  (n)  3w. 'w,  (n-1)  3w. • 3w J w,  (n- 1)  nj 

IK  IK  1  J  K 


(27) 
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32$ 


3w .  3w . 


i  j  k(n-l) 


Aw. 

J  (n) 


3<J> 


3w . 


(28) 


i  k(n-l) 


Eq.  (28)  is  used  for  the  Newton-Raphson  iterative  technique,  in 
which  wk(n  ^  represent  the  input  values  obtained  from  the  previous  solu¬ 
tion,  and  Aw_.  (n)  represent  the  correction  values  obtained  from  the  present 
iteration.  In  order  to  avoid  the  possible  divergent  solutions,  a  decelera¬ 
tion  coefficient  n  is  usually  assigned  such  that 


Wk(n)  "  Wk(n-1)  +  nAwk(n) 


(29) 


The  value  of  n  is  normally  restricted  to  0  <  n  £  1;  it  may  be  selected  as 

large  as  possible  if  the  functional  is  well  behaved.  Thus,  once  an  initial 

guess  wwm  is  provided,  the  iterative  relations  of  Eqs.  (28)  and  (29)  are 

then  used  to  obtain  an  improved  solution;  this  procedure  is  performed 

repeatedly  until  the  desired  accuracy  is  achieved.  Note  that  some 

32<J> 

restrictions  for  the  tangent  matrix,  [ - ],  are  required  in  order  to 

3w .  3w . 
i  J 

assure  the  existence  and  uniqueness  of  the  solution.  Further  discussions 
on  this  aspect  can  be  found  in  some  references  (e.g.,  [35]). 


2.  CONVERGENCE  REQUIREMENT 

In  order  to  ensure  the  convergence  of  the  obtained  solutions  during 
iterative  processes,  the  solution  must  be  checked  at  every  iteration  to 
satisfy  certain  requirements  according  to  a  convergence  criterion.  Since 
the  new  solution  is  controlled  by  the  size  of  the  deceleration  coefficient 
q  (Eq.  (29)),  the  provided  criterion  should  also  be  capable  of  determining 
an  optimal  value  of  the  deceleration  coefficient  in  order  to  have  fast 
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convergence  and,  hence,  reduce  the  computing  time.  Unfortunately,  no 
such  criterion  is  available  at  the  present  time,  although  efforts  have 
been  made  to  establish  one  in  the  past  years.  Instead,  the  following 
convergence  criteria  are  often  used  in  measuring  the  convergence  of  solu¬ 
tions. 

In  the  previous  studies  [17]- [21],  [29],  the  convergence  of  the  solu¬ 
tion  has  been  measured  by  the  quantity  l|  Aw ll/llw  II,  where  the  Euclidean  vector 
norm  is  defined  by 

/"Nxp  , 

llwll  =  /  l  w:  ,  (30) 

i-1  1 

where  N  =  total  number  of  nodal  points  and  p  =  number  of  unknowns  to  each 

node.  The  deceleration  coefficient  is  assigned  to  be  a  constant  and  the 

convergence  criterion  requires  that  the  error  norm  at  the  n-th  iteration 

(II Aw.  .ll/llw,  ,  J|)  be  less  than  that  at  the  previous  iteration.  This 
~(n)  ~(n-l) 

convergence  criterion,  with  a  suitable  selection  of  the  deceleration  co¬ 
efficient,  has  worked  out  well  for  the  solutions  in  previous  studies. 
However,  it  was  found  that  a  divergent  solution  may  be  indicated  according 
to  this  criterion  for  a  function  which  is  not  well  behaved,  although  the 
solution  is  converging  in  its  true  sense.  In  this  case,  the  criterion 
resulted  in  the  use  of  unnecessarily  small  values  of  the  deceleration 
coefficient  and  thus  increased  computing  time. 

Instead  of  the  measurement  defined  by  (30),  the  quantity  g  defined  by 
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I 

can  also  be  utilized  to  indicate  the  convergence  by  requiring  the  value 
of  g  evaluated  at  the  present  solution  (Eq.  (29))  be  less  than  that  of 
the  previous  one.  The  convergence  of  the  obtained  solutions  is  then 
measured  by  either  of  these  two  criteria  (Eqs.  (30)  and  (31)).  It  was 
found  that  a  combination  of  these  two  criteria  does  provide  an  adequate 
indication  of  convergence,  and  the  value  of  deceleration  coefficient  can 
also  be  selected  efficiently,  according  to  the  previous  information  of 
the  functional  behavior  and  the  convergence  requirement  at  the  current 
iteration. 

For  the  case  where  the  initial  guess  is  far  from  the  actual  one,  the 
convergence  requirement  by  II  Aw  ll/ll  w  II  is  very  helpful  in  preventing  the 
possible  diverging  solutions.  However,  once  the  obtained  iterative  solu¬ 
tion  is  very  close  to  the  actual  one  (say,  g  0.5),  the  measurement  by 
Eq.  (31)  and  its  convergence  requirement  can  provide  very  useful  informa¬ 
tion  in  selecting  an  efficient  value  of  the  deceleration  coefficient,  and, 
therefore,  provide  fast  convergence.  Note  that  situations  may  still 
happen  in  which  the  new  solution  obtained  does  not  satisfy  any  of  the 
criteria  above,  especially  when  a  very  bad  initial  guess  is  provided. 

This  difficulty  can  be  resolved  by  either  changing  the  initial  guess  or 
by  using  other  iterative  techniques  (for  example,  the  steepest  descent 
method  [36],  damped  Newton-Raphson  method  [33],  or  the  conjugate  gradient 
method  [37],  [38]).  Nevertheless,  it  was  found  that,  in  many  cases,  the 
bad  convergence  characteristics  are  often  due  to  the  improper  formulation 
itself  rather  than  the  choice  of  the  iterative  technique.  This  improper 
formulation  includes,  for  example,  the  undefined  value  of  the  frictional 
stress  (Eq.  (13))  for  the  neutral  region  along  the  die-workpiece  interface 
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for  the  ring  compression  case,  and  the  undefined  value  of  the  deviatoric 
stresses  in  Eq.  (3)  with  near-zero  effective  strain  rate  when  the  control 
volume  contains  rigid  portions.  The  necessary  modifications  for  the 
neutral  region  problem  will  be  discussed  later  in  Section  V.  The  method 
of  rigid  treatment  has  been  proposed  in  Section  II. 2.  Both  new  formula¬ 
tions  were  employed  in  the  present  investigations  and  the  convergent  solu¬ 
tions  were  therefore  obtained  without  difficulty. 

3.  DISCRETIZATION  FORMULATIONS  FOR  AXISYMMETRIC  AND  PLANE-STRAIN  CASES 

In  the  present  investigations,  the  nonlinear  stationary  value  of  the 
functional  4> ^  (Eq.  (14))  for  the  Lagrange  multiplier  and  the  minimization 
of  the  functional  <f>?  (Eq.  (15));  or  <t>3>  (Eq.  (18))  for  the  penalty 
constraint  were  solved  by  the  finite-element  method  and  the  iterative 
technique,  as  described  in  the  two  sections  above.  The  procedures  start 
with  the  approximation  of  the  functionals  above  by  Eq.  (33).  Let 
TT|ni'* ,  tt^,  be  the  integrals  defined  by  (32),  which  correspond 

to  distortional  energy,  dilatational  energy  of  the  Lagrange  form,  the 
penalty  form,  and  input  energy,  respectively.  All  these  integrals  are 
evaluated  at  typical  element  m.  These  integrals  are  expressed  by 


r(m)  . 
fl  '  J 


CTC  dV, 


,,0) 


.Cm)  _ 


,(m) 


A.c  dV 

V 


x(m)_(m).  ^(m)  _ 


,(m) 


£v  dV, 


(32) 


r(m)  _ 


Vm) 


e  dV. 
2  v 


r(n0  _  f 

F  -  J 


F*u  dS. 


(m) 
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where  denotes  the  volume  of  the  m-th  element  and  S^mj  denotes  the 

corresponding  surface  where  the  surface  traction  is  prescribed.  Note  that 
the  mean  stress  A  in  Eq.  (32)  is  assumed  to  be  uniform  (A  =  A^)  over  the 
entire  element.  The  functionals  and  <p 3  are  then  approximated  by 

M 

p,  =  J  +  A^tt^  -  (Lagrange  multiplier 

m=l  method) 

->  =  7  [tt^  +  1T-n1'*  -  ] .  (Penaltv  constraint) 

-  M= 1  1  J  h 

~  F  +  — ^  •  (Modified  penalty 

m=l  2V^‘  ]  constraint) 

where  M  is  the  total  number  of  elements. 

For  the  functional  the  independent  variables  are  the  nodal  point 
velocities  u^  (see  Eq.  (41))  and  the  element  mean  stress  A^;  the  system 
of  equations  to  be  solved  at  the  n-th  iteration  for  the  Newton- Raphson 
iterative  method  is  then  given  by 
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In  Eq.  (34)  the  tangent  matrix  is  banded  and  symmetric;  the  components  of 
this  matrix,  as  well  as  the  left-hand  side  vector  of  Eq.  (34),  involve 
only  the  nodal  point  velocities.  Hence,  only  the  initial  guess  of  the 


velocity  field  is  required  to  start  the  first  iteration;  the  value 

of  the  Lagrange  multiplier  (i.e.,  mean  stress)  can  be  obtained  from  the 
results  of  iterations. 

For  the  penalty  constraint  of  the  functional  <p^,  because  the  unknowns 
are  nodal  point  velocities  only,  the  system  of  equations  to  be  solved  for 
the  n-th  iteration  is 


M  m 

3  (, {■)  ♦  -  *<■))}  .Aufir  .  -  -  l  {^-(tt} 

m-1  3ua9u6  1  4  (n-1)  6(n)  m-1  oUa  1 


y  {, 


9  /_(m)  ,  _ (m)  (ra) 


♦  TT 1 


-  ) 


(35) 


The  expression  for  the  functional  <J>,  is  similar  to  Eq.  (34),  and  is  there- 

M 

fore  omitted.  Note  that  the  summations  (  J  )  in  Eqs.  (34)  and  (35)  denote 

m=l 

the  proper  assembling  procedure  by  adding  together  all  components  corres¬ 
ponding  to  the  same  unknown. 

The  following  discretization  formulations  are  made  at  the  elemental 
level  for  axisymmetric  cases,  in  which  a  four-nodal  quadrilateral  element 
is  considered.  The  necessary  modifications  of  these  formulas  for  plane- 
strain  cases  will  be  discussed  later.  The  actual  coordinate  system  (r,  z) 
for  a  typical  element  is  transformed  into  a  natural  coordinate  system 
(s,  t),  as  shown  in  Fig.  1.  This  natural  coordinate  system  is  defined  in 
a  way  such  that  s  and  t  vary  from  -1  to  1  within  each  element.  An  arbi¬ 
trary  point  (r,  z)  inside  the  element  can  therefore  be  expressed  in  terms 
of  the  natural  coordinate  (s,  t)  as  follows: 
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Fig.  1  Quadrilateral  element  and  natural  coordinate 
system. 


27 


(36) 


4 

r (s , t)  =  T  q  r  ; 

0-1  a  a 

4 

:(S>t)  =  ?  Vo  ’  (37) 

a=l  a  a  ' 

where  r  ,  z (a  =  1,...,4)  are  the  positions  of  the  four  surrounding  nodal 
points  of  a  typical  element  in  the  actual  coordinate  system;  and 

(1  *  s  s)  (1  +  t  t) 

=  - j - — —  (no  sum)  (38) 


are  called  the  transformation  functions.  The  values  of  s  ,  tQ  are 
tabulated  and  given  in  Table  1. 


Table  1 

Corresponding  values  of  s  and  t 
r  ®  a  a 

for  the  coordinate  transformation 


^"1  sT“| 

t 

a 

1  -1 

-i 

2  +1 

-i 

3  |  +1 

+i 

A  !  -1 

+i 

With  the  assumption  of  bilinear  distributed  velocity  fields,  the 

radial  and  axial  velocity  components  ur>  u^  at  any  arbitrary  point  inside 

the  element  can  then  be  determined  in  terms  of  the  surrounding  nodal 

4 

velocities.  This  is  expressed  by  u  =  £  N^u^  (Eq.  (22)),  or,  with  suitable 

a=l 

arrangement  of  the  representation,  expressed  by 


u. 

l 


8 

-  I 

a=l 


N .  u 
ra  a 


(i  =  1,2), 


(39) 
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(40) 


where  the  vector  {u^}  and  {u^}  are  defined  by,  respectively. 


and 


{u.}  =  {ur,u,}, 


T 

(u  }  =  {u  ,U,  ,U  ,U,  ,u  ,u„  ,U  ,U  }; 

1  "1  :  2  r5  “5  r4  “4 


(41) 


and  the  shape  function  [N^  ]  is  defined  by 


[N. 
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|V 

o, 

q2’ 

°,  q3. 

o, 

V 

0  1 

1 

Ip, 

V 

o. 

q2>  °» 

q5’ 

o. 

1 

V 

(42) 


The  Greek  index  a  denotes  the  appropriate  surrounding  nodal  point  vari¬ 
ables  at  the  elemental  level,  as  mentioned  before. 

The  stress  and  strain  components  in  Eqs.  (3)  and  (8)  can  be  repre¬ 
sented  by  the  proper  vector  forms,  such  that 


=  /  -r  e-g  , 


ana 


3  .  /  £  - , . 

*  -  '  7  u 


c '  , 


(43) 


(44) 


where  the  corresponding  values  of  the  vector  representations  are  given  by 

(45) 


and 


e  • 
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a! 
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A.  ,u.  , 
ij  J 


3  a 
7  —  £ 


i  ’ 


(46) 


where  [A^ ^ ]  is  the  appropriate  matrix  representation  of  the  differential 
operator  acting  on  the  velocity  component  {u^}.  For  axisymmetric  cases, 


{uJT  =  {ur>rz},  the  components  of  {e.},  {a!},  and  [A^]  are  defined  by 

{e±}T  =  {er.ez,e6./£r2}.  (47) 

(48) 


{a!}  =  (o', o', o'  Sl\ 


r’uz’uQ ’ ’ “c rzJ’ 


and 


[A.  .]  = 
1J 


_3_ 

3r 


J_  JL 

/2  9z 


3 

3: 


1  3 


/I 


3r 


(49) 


For  convenience,  the  appropriate  vector  representation  of  the 
Kronecker  delta  is  used  for  the  volumetric  strain  rate,  such  that 

ev  =  C-e,  (50) 

where  the  Kronecker  delta  C  is  defined  by 

(C}T  =  (1,1, 1,0}.  (51) 

By  substituting  the  expression  (40)  for  ur  and  u^  into  Eq.  (45),  the 
strain-rate  vector  is  then  expressible  in  terms  of  the  surrounding  nodal 
point  velocities  and  is  given  by 


=  B .  u 

io.  a 


and 


B. 

io 


A.  .N. 
11  jo 


(52) 


(53) 
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is  the  corresponding  matrix  representation  of  the  differential  operator 
A  acting  on  the  shape  function  N.  Results  of  these  differentiations  are 
obtained,  with  the  help  of  the  chain  rule,  by  converting  the  differentia¬ 
tions  from  the  (r,  z)  coordinate  to  the  (s,  t)  coordinate  according  to  the 
following  relationship: 


f 3 } 

r  3z 

rjL) 

l3r[ 

1  !  3t 

3s  j I 3s 1 

1  3 

'  W|  3r 

3rj]  3  [ 

L’3t 

3s-! 

where  J  is  the  determinant  of  the  Jacobian  of  the  coordinate  transforma¬ 
tion  and  is  expressed  by 

I J  i  -  Ir.li  .  ll.ll  (55) 

11  3s  3t  3t  3s  ' 


After  carrying  out  the  detailed  calculation,  the  values  of  the  matrix  B 
are  obtained  as  follows: 
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(60) 


In  the  expressions  above,  q  and  r  are  defined  by  Eqs.  (36)  and  (38);  r 

&  Up 

and  are  defined  by 


ct(3 


=  r  -  rc 
a  £ 


and 
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a0  a  3 

where  (r  ,  z  )  is  the  position  of  the  nodal  point  in  the  (r,  z)  coordinate 
system. 

According  to  Eqs.  (52),  (53),  and  (60),  the  effective  strain  rate  and 
volumetric  strain  rate  can  be  expressed  in  terms  of  the  surrounding  nodal 
velocities  as  follows: 


K  ,u  u- 
a3  a  8 


=  Q  u 

a 


(61) 

(62) 
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where  and  K^g  are  defined  by 


Kc,6  ■  8iaBia 


(63) 


and 
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„  =  w.B. 

a  1  ia 


(64) 


Inserting  all  the  corresponding  values  in  Eq.  (32)  and  performing 


the  necessary  differentiations,  the  first  and  second  derivatives  of  tt 
Tr^m^,  and  tt^  are  obtained.  These  are  given  by 
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Note  that  if  a  near-rigid  element  is  found  during  iterations,  the 
proposed  rigid  treatment  (Section  II. 2)  is  used.  In  this  case,  the  first 
and  second  derivatives  of  the  distortional  energy,  ,  become 


(m) 
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The  results  (65)  and  (66),  together  with  the  discretization  of  it ^ 
from  the  next  section,  are  substituted  into  Eqs.  (34)  and  (35).  With  the 
proper  assembling  procedures,  the  system  of  equations  necessary  for  the 
Newton-Raphson  iterative  technique  is  then  obtained.  This  system  of 
equation  is  used  to  obtain  a  new  solution,  and  the  solution  procedures 
are  repeated  until  the  desired  accuracy  is  achieved. 

Note  that  all  the  formulations  given  above  are  still  valid  for  plane- 
strain  cases,  if  the  following  changes  are  made:  (i)  the  (r,  z)  coordinate 
is  replaced  by  (x,  y)  coordinates;  (ii)  the  value  (i  =  1,...,4)  in 

Eq.  (59)  is  set  to  be  zero.  It  is  also  worth  noting  that  the  corresponding 
changes  in  the  numerical  evaluation  of  the  volume  integral  (Eq.  (65))  must 
be  made  for  plane-strain  cases.  Further  discussion  on  numerical  integra¬ 
tions  can  be  found  in  Section  1 1 1. 5. 


4.  DISCRETIZATION  FOR  TRACTION  PRESCRIBED  BOUNDARIES  AND  CALCULATION 
OF  RESULTANT  NODAL  FORCES 

In  general,  the  discretization  formulations  for  traction  prescribed 
boundaries  are  very  difficult  to  make  if  the  surface  traction  is  not 
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linearly  (or  uniformly)  distributed  over  the  surface  connecting  two 

consecutive  nodal  points.  For  this  reason,  the  following  discretization 

formulations  are  made,  under  the  assumption  of  linearly  distributed 

surface  traction  between  two  consecutive  nodal  points,  for  an  inclined 

boundary.  These  formulations  were  used  for  the  cases  of  extrusion  and 

drawing,  in  which  the  interfacial  frictional  stress  is  assumed  as  a 
Y0 

constant  (f  =  m  — )  and  its  direction  is  known. 

/3 

The  inclined  boundary  is  shown  in  Fig.  2,  in  which  the  surface  makes 
an  angle  a  with  the  horizontal  axis.  The  positions  (r^  z^)  and  (r2,  z2) 
are  assigned  to  nodals  1  and  2,  respectively,  for  a  typical  element  along 
this  boundary.  Normal  and  shear  stresses,  (pj,  t.)  and  (p2>  t2)  are 
given  at  these  points.  With  the  assumption  of  linearly  distributed 
traction,  the  prescribed  stresses  in  the  r  and  z  directions,  at  an  arbi¬ 
trary  point  (r,  'z)  along  this  surface,  are  given  by 
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sin  a. 


For  the  quadrilateral  element  used  in  the  present  investigations, 
the  velocities  have  bilinear  distributions  inside  the  element  and  are 
linear  along  the  boundaries.  For  any  arbitrary  point  along  the  surface, 
the  velocities  (ur,  uz)  are  expressible  in  terms  of  the  nodal  velocities 


u  tr,  -  r)  u  (r,  -  r) 

1  1  j  1 
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u  = 
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Inserting  Eqs.  (67a)  and  (67b)  into  Eq.  (32d)  and  performing  the 
necessary  integration,  the  following  results  are  then  obtained  for  the 
integral  for  plane  cases  and  axisymmetric  cases,  respectively: 
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i axisvmmet ric  cases: 


(69) 


where  W  is  the  width  of  the  workpiece  in  plane  cases,  and  Ar,  Az  are  de¬ 
fined  by  Ar  =  r2  -  ^  and  Az  =  z2  -  z^,  respectively. 
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In  view  of  Eqs.  (68)  and  (69),  the  first  derivatives  of  which 

r 

are  equivalent  to  the  corresponding  nodal  forces,  are  all  constants; 
and  the  second  derivatives  are  all  zeros.  These  results  are  substituted 
into  Eqs.  (34)  and  (35)  for  the  formulation  of  Newton-Raphson  iterations. 

There  is  another  type  of  boundary  condition  in  which  the  frictional 
stress  depends  not  only  on  the  direction,  but  also  on  the  magnitude  of 
the  velocity  of  the  deforming  body  relative  to  the  die,  such  as  the  neutral 
flow  problems  described  in  Section  V.  The  discretization  procedures  are 
very  complicated  and  the  formulations  require  very  tedious  manipulations. 
Detailed  discretization  formulations  for  this  kind  of  boundary  conditions 
are  given  in  Section  V. 

For  nodal  points  along  the  velocity  prescribed  boundary,  the  nodal 

velocity  u^  is  equal  to  the  prescribed  value  and  the  obtained  value  of 

Au^  must  always  be  equal  to  zero.  This  can  be  accomplished  by  modifying 

the  equation  corresponding  to  Au^  in  Eq.  (34)  or  (35),  such  that  this 

equation  becomes  Au^  =  0  during  every  iteration.  The  obtained  velocity 

field  and  the  value  of  the  Lagrange  multiplier  are  then  used  to  calculate 

the  resultant  nodal  force  corresponding  to  this  velocity  prescribed  nodal 

point,  according  to  ffl  =  £  [ir^  +  X  ]  for  the  Lagrange 

m=l  ua 


multiplier  constraint,  and  fQ  =  J  +  17 ^  ] 

m=l  °  ~ 


M 


or  \  (TT^m^)^])  for  the  penalty  function  constraint. 

m=l  dua  2vlm'  v 


5.  DISCUSSION  OF  THE  NUMERICAL  INTEGRATION 

Throughout  the  present  investigation,  the  Gaussian  quadrature  formulas 
are  used  for  volume  integration.  For  a  given  function,  H(r,  z)  (or  H(x,  y) 
in  plane  cases)  the  integration  is  given  by 
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(plane  cases) 
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where  W  is  the  width  of  the  workpiece  in  plane  cases  and  |j|  is  the 
determinant  of  the  Jacobian  of  the  coordinate  transformation. 

The  2x2  Gaussian  quadrature  gives  the  following  evaluation  of  these 
volume  integrations: 


2tt  }  P±  [H  *r  •  J  ]  j 
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H  dV  = 


s=s. ;  t=t . 
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i  =  l 


s=s . ;  t=t . 

l  l 


where  is  the  weighting  factor,  which  is  equal  to  1  for  all  four  inte¬ 
gration  points.  The  locations  of  these  integration  points  are  tabulated 
and  given  in  Table  2. 

For  an  integral  with  its  integrand  expressible  by  a  polynomial  of  s 

and  t,  the  exact  evaluation  of  this  integral  can  be  made  by  the  2x2 

3  3 

Gaussian  quadrature  for  orders  up  to  s  and  t  .  In  the  present  studies, 
this  2x2  numerical  integration  is  applied  to  all  volume  integrals  except 
7r^m)  (tt ^  z"  dV) .  Note  that  a  reduced  lxi  integration  must  be 

3  3  Jv(m)  1  V 

used  for  ir^11^  in  order  to  avoid  the  possibility  of  over-constraint  due  to 
the  2x2  numerical  integrations.  The  reasons  are  clearly  shown  below. 
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Table  2 


Locations  of  integration  points 
for  2*2  Gaussian  quadrature  formulas 
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Consider  a  typical  element  V^;  the  volumetric  strain  rate 
inside  this  element  is  given  by 


e  = 
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where  a^,  b^,  c^  are  all  independent  of  (s,  t) ;  and  the  detailed  expres¬ 
sions  of  these  constants  can  be  found  from  Eqs.  (62)  and  (56),  which  are 
omitted  here.  For  the  cases  of  plane-strain  and  plane-stress,  the  exact 
evaluation  of  the  integral  corresponding  to  dilatation  energy  in  the 
Lagrange  multiplier  formulation  is  given  by 


4m)  =  A(m) 


.  ,  ,  8 
£  dV  =  4A  m  (  y 
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a  u  )  , 
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(73) 


which  can  also  be  obtained  from  either  a  2x2  numerical  integration  or  a 
reduced  1*1  integration.  Note  that  in  the  lxl  integration,  the  integration 
point  is  located  at  the  centroid  of  the  element  (i.e.,  s  =  0,  t  =  0)  and 
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the  weighting  factor  is  4.  Now  consider  the  penalty  function  constraint 
tt^  ;  the  integral  can  be  written  in  terms  of  the  local  coordinates  s  and  t: 


-1  J-l  2  (aa  +  V  +  V^ct^  (a3  +  b6s  +  ds  dt‘ 


With  reduced  lxl  numerical  integration,  the  result  of  the  integral  above 
is  given  by 
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where  A  is  the  area  of  the  element.  Clearly,  the  derivatives  of  the  equa¬ 
tion  above  (75)  is  equivalent  to  those  derivatives  of  Eq.  (73)  if  the 
mean  stress  A^  is  made  according  to  (17)  (i.e.,  A  ^  =  ££v)  with  the 

value  of  e  evaluated  at  the  centroid  of  the  element, 
v 

Hence,  when  the  minimum  solution  is  made,  the  1*1  integration  of  the 

penalty  constraint  tt^  gives  the  desired  incompressibility  condition, 

which  is  equivalent  to  that  of  the  Lagrange  multiplier  constraint.  This 

is,  however,  no  longer  true  if  the  2><2  numerical  integration  is  used  in 

f  C  •  2 

the  evaluation  of  the  integral  y  e  dV,  because  the  constants  b  and  c 

y  ^  v  a  a 

in  Eq.  (73)  and  the  reduced  one-point  integration  of  (74)  require  the  zero 
volumetric  strain  rate  be  satisfied  only  at  the  centroid  of  the  element; 
but  the  2x2  integration  requires  the  zero  values  of  ey  be  satisfied  at  all 
four  integration  points.  This  over-constraint  requirement  induced  by  the 
2x2  integration  will  give  an  incorrect  solution;  the  lxl  reduced  integra¬ 
tion  is,  therefore,  highly  recommended  for  the  evaluation  of  the  penalty 
term,  Tr^  . 
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It  is  to  be  noted  that  if  the  modified  penalty  function  — 7— 

2v(m)  v 

(Eq.  (33c)  is  used,  the  reduced  integration  is  no  longer  necessary  because 

•2 

the  quantity  is  not  presented  in  the  formulation.  Thus,  the  modified 
penalty  function  can  be  applied,  without  any  restriction,  to  any  problems, 
including  axisymmetric  and  three-dimensional  cases. 
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Section  IV 


MULTIPASS  EXTRUSION  AND  DRAWING 


1.  INTRODUCTION 

Many  bar  drawing  and  extrusion  processes  in  practice  are  multipass 
operations.  Internal  fracturing  (center  bursting)  in  these  processes 
usually  occurs  after  some  area  reduction  through  several  passes.  In  order 
to  predict  and  prevent  the  occurrence  of  such  defects,  it  is  necessary  to 
determine  the  deformation  characteristics  involved  in  multipass  operations 
The  influence  of  the  deformation  history  on  the  current  flow  pattern, 
among  others,  is  a  unique  feature  of  the  problem. 

Since  the  detailed  mechanics,  particularly  the  stress  and  strain 
paths  at  a  critical  site  of  a  deforming  material,  are  the  required  informa 
tion,  the  method  of  analysis  should  be  capable  of  accurately  determining 
not  only  overall  quantities  involved  in  metalworking  processes,  such  as 
forming  loads,  but  also  stress  and  strain  distributions  during  deformation 
Furthermore,  work-hardening  and  distributed  flow  stress  due  to  previous 
passes  must  be  taken  into  consideration  in  the  analysis.  The  finite- 
element  approach  comes  close  to  fulfilling  these  requirements. 

In  the  present  investigation  the  rigid-plastic  finite-element  method 
was  used  as  an  analysis  technique.  The  effect  of  work-hardening  was 
examined  in  single-pass  operations  and  the  effect  of  distributed  flow 
stress  on  the  flow  pattern  was  revealed  by  analyzing  two-pass  drawing  and 
extrusion. 
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2.  METHOD  OF  ANALYSIS 

The  steady-state,  multipass,  axisymmetric  extrusion,  and  drawing 
processes  were  analyzed  by  the  rigid-plastic  finite-element  method  with 
the  modified  penalty  function  constraint,  Eq.  (18).  Numerical  integra¬ 
tions  of  the  modified 

one-point  integration  evaluated  at  the  element  center.  Fig.  5  compares 
the  results  of  2*2  integration  and  1*1  integration  for  the  penalty  term 
for  the  case  of  extrusion  without  friction.  As  Fig.  5  shows,  the  reduced 
one-point  integration  resulted  in  a  smoother  distribution  of  the  stress 
along  the  axis  of  symmetry,  with  its  magnitude  approximately  equivalent 
to  the  mean  value  of  the  stress  obtained  by  2x2  integration. 

Integration  of  strain  rates  were  made  according  to  the  scheme  dis¬ 
cussed  in  the  previous  AFML  report  [40],  For  work-hardening  materials  at 
the  first  pass,  the  effective  strain  was  assumed  zero  at  the  starting  point 
of  the  flow  line.  For  the  second  pass  (or  any  pass  other  than  the  first 
pass),  the  resulting  effective  strain  of  the  previous  pass  was  used  as 
the  initial  value  of  the  effective  strain  at  the  starting  point  of  the 
flow  line.  Effective  strains  were  calculated  at  every  iteration,  and  the 
iteration  was  performed  repeatedly  until  the  desired  accuracy  was  achieved. 
Note  that  the  stress  distributions  obtained  from  computations  differ  from 
the  actual  distribution  by  a  uniform  hydrostatic  component.  Accordingly, 
the  necessary  corrections  were  made,  using  the  technique  described  in  [40] . 

3.  COMPUTATIONAL  CONDITIONS 

Figure  3  shows  the  mesh  systems  and  the  boundary  conditions  used  for 
the  analysis.  Both  systems  consist  of  135  elements  interconnected  at  160 
nodals.  The  process  conditions  for  extrusion  and  drawing  under  which  the 


penalty  term. 


2V 


(m) 


,(m) 


e  dV]  ,  were  made  by 
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computation  was  carried  out  are  summarized  in  Tables  3  and  4,  respectively. 


Table  3 

Extrusion  process  conditions  for  computation 
Die  semi-cone  angle:  a  =  45° 
Reduction  ratio:  R^/R^  =  1-6 

Die  workpiece  interface 


frictional  stress  f 


Material 


1st  pass 


Non-work- hardening 
O  =  Y0 

Pure  aluminum 


i  0 

I 

!  0.4Y 

i 


0 


-  v  , ,  £ _  .0.54 

v  =  Y0'1  +  0.519 


YQ  =  106 . 86MN/m" 


I  'J-4Vo 


2nd  pass 


0 

0.4Y 


0 


Table  4 


Drawing  process  conditions  for  computation 


Die  semi-cone  angle:  a  =  6° 
Reduction  ratio:  R^/R^  =  1.1 


Die  workpiece 

interface 

frictional  stress  f 

Material 

1st  pass 

2nd  pass 

Non- work -hardening 

5  -  Y0 

°-4Yo 

- 

Aluminum  alloy 

2024-T351 

0=Y(1+  °  0 . 1 71 

V  0.0102-* 

°-4Y0 

°-4Yo 

Y0  =  309.4MN/m 

45 


and  the  boundary  conditions 
ation  for  (a)  extrusion 


AL  alloy  2024  1 351 


Fig.  4  Stress-strain  curves  of  the  two 
materials  used  in  computation. 


Fig.  5  Comparison  of  one-point  and  2x2 

Gaussian  integration  of  the  penalty 
term  (axial  stress  vs.  total  strain 
along  the  axis  of  extrusion; 
material:  non-work -hardening) . 
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Stress-strain  curves  for  pure  aluminum  and  aluminum  alloy  2024-T351  are 
shown  in  Fig.  4.  The  convergence  is  measured  by  the  norm  defined  by 


j|Aui|  /  £  2  -,2i 

-  — ,  where  |iui|  =  /  l  Uu  )“  +  (u.  )  ) 

'  i  =  l  i  “i 

N  being  the  total  number  of  the  nodal  points.  The  iteration  continues 
for  each  solution  until  the  accuracy  of  Au  /  u  <_  0.00005  is  reached. 
The  number  of  iterations  needed  to  achieve  this  convergence  depends  on 
the  initial  guess  of  the  velocity  field.  The  average  number  required 
for  the  final  solution  was  20  30  iterations. 

4.  EXTRUSION 

The  average  extrusion  pressures  computed  for  all  of  the  cases  given 
in  Table  3  are  summarized  in  Table  5. 

Table  5 


Average  extrusion  pressures 

Die  semi-cone  angle:  a  =  45° 

Reduction  ratio:  R./R„  =1.6 
i  0 


Material 

!  f  =  o 

!  f  =  °-4V0 

Non- work- hardening 

iW'o  • 

1 — 

1  1.823 

Pure  aluminum 

| 

1st  pass 

2.218 

2.642 

2nd  pass 

— -  j  —  —B  — — —  -  .  i  .  .  - 

2.982  | 

3.425 

We  shall  now  discuss  the  detailed  solutions  in  terms  of  velocity  distribu¬ 
tions,  grid  distortions,  and  stress  and  strain  distributions. 


The  general  trend  of  the  velocity  distributions  is  the  same  in  all 
cases  and  is  in  agreement  with  the  measurements  in  the  visioplasticity 
study  of  axisymmetric  extrusion  of  lead  [41].  In  Fig.  6  the  effects  of 
work-hardening  are  shown  by  comparing  the  velocity  distributions  for  a 
non-work-hardening  material  and  for  commercially  pure  aluminum  under  the 
interface  friction  of  f  =  0.4Yq.  The  two  velocity  components  are  plotted 
as  functions  of  the  radial  coordinate  at  various  locations  in  the  axial 
direction.  The  distribution  patterns  are  the  same  for  both  cases  but  the 
magnitudes  differ.  The  velocity  distributions  during  the  first  pass  and 
the  second  pass  in  frictionless  extrusion  are  compared  in  Fig.  7.  A  close 
examination  of  Figs.  6  and  7  shows  that  the  distributions  during  the 
second  pass  are  closer  to  those  for  a  non-work-hardening  material. 

Detailed  differences  of  deformation  characteristics,  due  to  material 
properties  and  multipass  operations,  are  more  clearly  revealed,  for  example, 
in  grid  distortions.  The  steady-state  grid  distortion  patterns  are  com¬ 
pared  for  non-work-hardening  and  work-hardening  cases  in  Fig.  8  and  for 
the  two  passes  in  Fig.  9.  It  is  seen  in  Fig.  10  that  distortion  of  grid 
lines  that  are  originally  perpendicular  to  the  axis  of  the  workpiece  is 
less  for  non-work-hardening  materials  than  that  for  work-hardening  materials. 
A  comparison  of  grid  distortions  for  two  passes,  given  in  Fig.  9,  shows  that 
distortion  of  grid  lines  for  the  second  pass  is  less  than  for  the  first 
pass,  confirming  the  observation  in  velocity  distributions  that  the  deforma¬ 
tion  pattern  during  the  second  pass  is  closer  to  that  for  non-work -hardening 
materials.  An  important  feature  of  grid  distortions  can  be  seen  from  the 
vertical  grid  distortion  at  the  exit,  as  plotted  in  Fig.  10.  For  friction 
extrusion,  grid  distortion  is  greatest  for  the  first  pass  of  work-hardening 
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Fig.  6  Distributions  of  velocity  components 
a  =  45°,  R./R0  =1.6,  f  =  0.4Y0; 

materials:  non-work-hardening  ( - ) 

and  pure  aluminum,  1st  pass  ( - ). 


Fig.  7  Distribution  of  velocity  components 
a  =  45°,  Ri/RQ  =  1.6,  f  =  0; 

material:  pure  aluminum,  1st  pass  (■ 
2nd  pass  ( - ) . 


PURE  ALUMINUM 


NON-WORKHARDENING 


Fig.  8 


Grid  distortion  for  non¬ 
work  -hardening  materials 
and  pure  aluminum  (f  =  0.4Y 


0' 


Fig.  9  Comparison  of  grid  distortion 
in  pure  aluminum  for  1st  and 
2nd  passes  (f  =  0)  . 
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materials,  less  for  the  second  pass  of  work-hardening  materials,  and 
least  for  non-work-hardening  materials.  This  order  is  different,  however, 
when  interface  friction  is  present,  because  the  magnitude  of  the  frictional 
stress  relative  to  the  flow  stress  of  incoming  material  is  less  for  the 
second  pass  of  work-hardening  materials,  in  comparison  to  other  cases. 

This  is,  in  effect,  a  reduction  of  friction,  thereby  reducing  the  total 
grid  distortion. 

The  effective  strain  rates  were  integrated  along  the  flow  lines  to 
yield  the  effective  strain  distributions  and  are  shown  in  Fig.  11.  The 
figure  also  shows  the  boundaries  between  the  rigid  and  plastically  deform¬ 
ing  regions  by  a  near- zero  strain-rate  contour.  The  rigid-plastic  boundaries 
at  the  exit  are  almost  identical  in  shape  and  location  for  all  cases.  How¬ 
ever,  the  size  of  the  deformation  zone  is  larger  for  work-hardening  mater¬ 
ials  at  the  first  pass  but  reduces  at  the  second  pass  and  becomes  almost 
the  same  as  that  for  non-work-hardening  materials.  It  is  of  interest  to 
note  that  the  distribution  of  additional  strains  during  the  second  pass 
is  very  close  to  that  for  non-work-hardening  materials.  So  far,  all  indi¬ 
cations  with  regard  to  the  effect  of  distributed  flow  stress  on  the  incom¬ 
ing  material  are  that  the  deformation  characteristics  during  the  second 
pass  are  different  from  those  during  the  first  pass  due  to  the  non-uniform 
flow  stress  distributions  in  the  incoming  material  and  are  closer  to  those 
for  non-work-hardening  materials.  The  reason  for  this  is  not  obvious. 
Examining  the  strain  distributions  and  corresponding  flow  stress  distribu¬ 
tions  at  the  exit  section  shown  in  Fig.  12,  we  can  see  that  non-uniformity 
of  flow  stress  distribution  after  the  first  pass  is  appreciable  due  to 
the  effect  of  work-hardening,  but  that  the  increase  in  flow  stress  during 
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Fig.  11  The  effective  strain  distribu¬ 
tions  in  extrusion. 
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t»t  PASS 


Fig.  12  The  effective  stress  and  the 

effective  strain  distributions 
across  a  section  of  an  extruded 
bar. 
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the  second  pass  is  almost  constant  across  the  section.  During  the  second 
pass,  the  two  effects  are  combined.  These  are  the  effects  of  further 
work-hardening  and  the  effect  of  initial  nonuniform  flow  stress  distribu¬ 
tions.  Due  to  work-hardening  during  the  first  pass,  the  manner  in  which 
the  material  work-hardens  during  the  second  pass,  however,  is  different 
because  the  state-of-material  points  are  located  at  different  points  on 
the  stress-strain  curve.  These  combined  effects  appear  to  be  influential 
factors  in  determining  the  deformation  characteristics  in  multipass  extrusion. 

The  distribution  of  hydrostatic  pressure  (-°m/^)  were  examined  in  Fig. 

13  for  frictionless  extrusion.  It  was  found  that  the  distribution  patterns 
are  the  same  for  all  cases  shown  except  for  some  differences  in  magnitude. 

It  was  observed  that  along  the  die  workpiece  interface,  the  hydrostatic 
pressure  shows  a  minimum  at  some  distance  from  the  exit.  This  feature  is 
seen  also  in  the  normal  pressure  distributions  along  the  die  interface 
shown  in  Fig.  14.  For  non-work-hardening  materials,  the  pressure  is  high¬ 
est  at  the  entrance  and  it  decreases  toward  the  exit,  increasing  again 
near  the  exit.  For  work-hardening  materials,  the  trend  is  the  same  but 
the  pressure  at  the  exit  is  slightly  higher  than  that  at  the  entrance. 

For  the  second  pass,  the  magnitude  increases  again  with  the  pressure  great¬ 
est  at  both  the  exit  and  the  entrance. 

The  critical  site  for  the  occurrence  of  center  bursting  in  extrusion 
is  along  the  s  of  extrusion.  The  stress  and  strain  histories  are  im¬ 
portant  for  appl  ng  a  ductile  fracture  criterion.  Fig.  15  shows  these 
stress  and  strain  histories  along  the  axis  of  extrusion  for  the  friction¬ 
less  case.  The  axial  stress  becomes  tensile  over  a  certain  range  of  de¬ 
formation  and  a  larger  tension  zone  is  obtained  for  non-work-hardening 
materials  than  that  for  work-hardening  materials.  The  tensile  stress 
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increases  during  the  second  pass.  It  was  found  that  the  interface  fric¬ 
tion  has  the  effect  of  decreasing  the  axial  stress,  thus  suppressing  the 
possibility  of  producing  center  burst.  These  results  are  qualitatively 
in  agreement  with  experimental  observations.  Fig.  16  shows  the  variations 
of  axial  stress  along  the  center  axis  for  the  case  of  constant  interface 
friction  with  f  =  0.4Y^.  As  mentioned  earlier,  the  magnitude  of  the 
frictional  stress  (f  =  0.4Yq)  relative  to  the  flow  stress  (5)  is  less  for 
the  second  pass  than  that  for  the  first  pass.  Therefore,  due  to  the  effects 
of  this  reducing  frictional  stress,  the  amount  of  tensile  stress  increased 
during  the  second  pass  is  greater  for  the  friction  case  (Fig.  16)  than  that 
for  the  frictionless  case  (Fig.  15). 


5.  DRAWING 

The  drawing  process  is  characterized  by  smaller  die  angles  and  smaller 
area  reduction,  when  compared  with  those  for  extrusion.  The  computed  aver¬ 
age  drawing  stresses  are  listed  in  Table  6. 

Table  6 

Average  drawing  stresses 
Die  semi-cone  angle:  a  =  6° 


Reduction  ratio: 

VR0  -1.1 

Material  j  f  =  Q.4Yq 

i 

Non-work-hardening  i 

a  A'  =  0.906 
ave  0 

Aluminum  allov  2024-T351  i 

1st  pass 

1.042 

2nd  pass 

1. 118 
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15  Stress  and  strain  history  of  a  material  point  along 
the  axis  of  extrusion  (f  =  0). 


lit  PASS  ,2nd  PASS 


Grid  distortions  at  the  exit  are  given  in  Fig.  17.  In  comparison 
with  Fig.  10  for  extrusion.  Fig.  17  shows  that  the  grid  lines  for  drawing 
form  the  double-peak  distortion  in  the  drawn  part  for  non-work-hardening 
materials.  For  work-hardening  materials,  this  tendency  lessens  for  the 
first  pass,  but  increases  again  for  the  second  pass,  suggesting  that  the 
deformation  characteristics  during  the  second  pass  become  similar  to  that 
of  non-work-hardening  materials. 

The  strain  distributions  (Fig.  18)  for  the  three  cases  are  similar 
to  those  for  extrusion,  as  shown  in  Fig.  11.  The  variation  of  effective 
strain  distributions  in  drawn  bars  (Fig.  19)  is  also  similar  to  that  in 
Fig.  12  for  extrusion  except  that  the  variation  is  comparatively  small 
due  to  a  small  die  angle  and  small  reductions.  In  drawing  it  is  observed 
also  that  the  increase  in  flow  stress  during  the  second  pass  is  almost 
uniform  across  the  bar  section.  The  hydrostatic  component  distributions 
are  given  in  Fig.  20.  Also,  the  rigid-deforming  boundaries  are  shown  in 
the  figure.  The  deformation  zone  size  increases  in  the  first  pass  of 
work -hardening  materials,  compared  to  that  of  non-work-hardening  materials, 
but  decreases  in  the  second  pass  and  becomes  almost  identical  to  that  for 
non-work -hardening  materials.  It  is  also  concluded  that  the  strain  distribu¬ 
tions  during  the  second  pass  are  closer  to  those  for  a  non-work-hardening 
material,  although  the  difference  between  those  for  the  first  pass  of 
work-hardening  materials  and  of  non-work-hardening  materials  are  not  too 
great.  The  hydrostatic  component  distributions  revealed  that  the  region 
under  the  die  is  in  the  compressive  hydrostatic  stress  state,  while  the 
core  and  exit  regions  are  in  the  tensile  hydrostatic  stress  state.  The 
stress  distributions  are  similar  with  some  differences  in  magnitude.  The 
die  pressure  distributions  are  greatest  at  the  entrance  and  minimum  at  the 
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_ _ _ NON-WORK  HARDENING 
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17  Grid  distortion  at  exit  section 
in  drawing  (f  =  0.4YQ). 
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Fig.  18  The  effective  strain  distribu¬ 
tions  in  drawing. 


Fig.  20  Hydrostatic  pressure  distribu 
tion  in  drawing. 


halfway  point  of  the  interface  for  all  three  cases,  as  shown  in  Fig.  21. 

The  stress  and  strain  histories  along  the  axis  of  drawing  are  plotted 
in  Fig.  22.  The  axial  stress  is  greater  for  non-work-hardening  materials 
than  for  work-hardening  materials  during  the  first  pass,  but  increases 
during  the  second  pass.  The  peak  stress  appears  near  the  entrance,  while 
this  peak  is  observed  near  the  exit  for  extrusion. 

6.  CONCLUDING  REMARKS 

The  rigid-plastic  finite-element  method,  with  the  penalty  function 
constraint,  has  been  successfully  applied  to  the  analysis  of  multipass 
bar  drawing  and  extrusion.  The  velocity  distributions,  grid  distortions, 
die  pressure,  and  stress  and  strain  distributions  were  obtained  in  detail 
for  two-pass  extrusion  and  drawing  of  a  work-hardening  material.  Solutions 
were  also  obtained  for  non-work-hardening  materials.  Deformation  char¬ 
acteristics  during  the  second  pass  were  generally  closer  to  those  of  non¬ 
work-hardening  materials,  in  comparison  with  first-pass  deformation,  for 
both  extrusion  and  drawing.  The  axial  stress  along  the  axis  of  extrusion 
and  drawing  was  found  to  be  greater  in  tension  for  the  second  pass  than 
that  during  the  first  pass,  suggesting  an  increasing  possibility  of 
center  burst  formation  in  multipass  processes. 
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Section  V 


RING  COMPRESSION 


1.  INTRODUCTION 

When  a  short  ring  specimen  is  compressed  between  two  flat,  parallel 
platens  (Fig.  23(a)),  the  diameter  of  the  inner  surface  either  increases 
or  decreases  as  the  height  of  the  specimen  is  reduced,  depending  on  the 
friction  condition  at  the  interface.  The  inner  diameter  of  the  ring 
is  increased  if  the  interface  friction  is  low  (Fig.  23(b))  and  it  is 
decreased  if  this  friction  is  high  (Fig.  23(c)).  Because  the  change 
in  internal  diameter  of  the  compressed  ring  is  sensitive  to  friction 
at  the  die-workpiece  interface,  ring  compression  has  been  widely  used 
as  a  test  to  evaluate  the  friction  condition  in  metal-working  processes. 

In  order  to  demonstrate  the  friction  condition  quantitatively,  how¬ 
ever,  the  relationship  between  the  geometrical  change  of  the  workpiece 
and  the  friction  condition  at  the  work  tool-workpiece  interface  must  be 
established.  This  aspect  of  ring  compression  has  been  the  subject  of 
many  investigations  in  the  past  [42]- [55].  Most  of  the  mathematical 
methods  are  based  on  an  approximate  stress  analysis  or  on  the  upper- 
bound  analysis  using  relatively  simple  velocity  fields.  Earlier  analy¬ 
sis  assumed  that  the  cylindrical  surface  of  the  ring  specimen  remains 
cylindrical  during  compression.  Explicit  expressions  of  neutral  point 
radius  and  corresponding  geometrical  changes  were  obtained  by  the 
stress  analysis  for  rigid-perfectly  plastic  materials  with  constant 
coefficient  friction  [42]  or  with  constant  interface  frictional  stress 
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[45].  Under  the  same  assumptions,  the  upper-bound  solutions  were 
given  for  constant  frictional  stress  [43],  [44],  DePierre  and  Male 
[47]  attempted  to  include  the  effect  of  work-hardening  by  assuming 
that  the  ratio  of  the  frictional  stress  to  the  current  flow  stress 
remains  constant  during  the  process.  They  observed  in  comparing  theory 
with  experiment  that  the  test  results  under  high  friction  conditions 
exceeded  the  theoretical  maximum  curve  where  the  frictional  stress  is 
equal  to  the  shear  yield  stress.  This  discrepancy  was  also  found  in 
a  subsequent  report  [48]  in  predicting  material  flow  stress  using  ring 
compression.  It  was  suggested  that  the  discrepancy  may  be  due  to 
severe  barreling  of  the  ring-free  surface  under  conditions  of  high  fric¬ 
tion. 

The  upper-bound  solutions,  including  bulge  formation,  were  intro¬ 
duced  by  Avitzur  [49],  assuming  exponential  cusp-shaped  bulge,  and  by 
Liu  [51],  for  cosine-shaped  bulge.  In  addition  to  the  bulge  formation, 
material  dependency  of  deformation  was  taken  into  the  analysis  by  Lee 
and  Altan  [52].  They  utilized  a  finite  mesh  system  for  computing  the 
flow  stress  distribution  and  minimized  the  total  energy  rate  according 
to  the  provided  admissible  velocity  field. 

While  these  theories  provided  useful  information  in  practice,  more 
accurate  and  realistic  solutions  are  needed  for  wider  and  more  meaning¬ 
ful  applications  of  the  ring  test  in  industry. 

2.  RIGID-PLASTIC  FINITE-ELEMENT  METHOD 

A  unique  feature  of  deformation  in  ring  compression  (also  in  rolling 
and  forging)  is  that  there  exists  a  point  (or  a  region)  along  the  die- 
workpiece  interface  where  the  velocity  of  the  deforming  material 


relative  to  the  die  velocity  becomes  zero  (neutral  flow  point) .  The 
frictional  stress  usually  changes  its  direction  at  the  neutral  flow 
point,  but  the  location  of  this  point  is  not  known  a  priori.  In 
implementing  the  finite-element  method  for  the  analysis  of  this  type 
of  problem,  difficulties  always  exist  in  the  formulation  of  boundary 
conditions  along  the  die-workpiece  interface.  Recently,  Matsumoto, 

Oh,  and  Kobayashi  [56]  reformulated  the  boundary  conditions  according 
to  the  direction  of  relative  velocities  at  the  interface,  and  solved 
the  problem  of  ring  compression  with  low  friction.  However,  a  diffi¬ 
culty  was  encountered  for  the  case  of  high  friction,  because  of  the 
development  of  the  neutral  flow  region  (rather  than  the  neutral  flow- 
point)  .  In  the  present  investigation,  a  variational  formulation  was 
developed  for  the  case  where  the  frictional  stress  is  dependent  on 
the  relative  velocity  at  the  die-workpiece  interface.  This  formula¬ 
tion,  when  interface  friction  is  low,  is  almost  identical  to  that  devel¬ 
oped  by  Matsumoto  et  al.,  and  permits  the  treatment  of  neutral  flow 
region  along  the  interface.  The  finite-element  solution  of  ring 
compression  was  then  obtained  under  various  friction  conditions  for 
several  materials.  Some  of  the  results  were  compared  with  the  pre¬ 
viously  published  theoretical  solutions  and  experimental  results. 

A .  Problem  formulation 

Consider  a  body  of  volume,  V,  with  the  traction,  F,  prescribed  on 

a  part,  Sp,  of  the  surface  and  the  velocity,  V,  prescribed  on  surface, 

S  .  Let  S  be  the  remainder  of  the  surface  where  the  frictional  stress, 
u  c 

f,  acts  and  the  die  velocity,  \jP  ,  is  given.  The  actual  stress  and 
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velocity  fields  satisfy  the  following  relations: 


(a)  Equilibrium  conditions: 


a.  .  .  =  0  (neglecting  the  body  force); 
iJ.J 


(76) 


(b)  Compatibility  conditions: 


e.  .  = 
ij 


=  4(v.  .  +  v.  . )  and  e  =  v. 

2  i,j  j , l  v  i,i 


=  0- 


(77) 


(c)  Stress-strain  rate  relationship: 


,  _  2  __0__  • 

Jij  -  3  |  cij  : 


(78) 


where  a!,  is  the  deviatoric  stress,  c  and  t  are  the 
ij 

effective  stress  and  the  effective  strain  rates,  respec¬ 


tively,  and  defined  by 
(d)  Boundary  conditions: 


.AT— 


r 


•  • 


—  v  -r  o :  .  o !  .  and  z  —  v  —  z .  -  c .  . 
-  11  l]  o  11  11 


a.  .n.  =  F  .  on  Sr. 

1J  1  1  F 


v.  =  V.  on  S  . 
11  u  ’ 


a.  . n.  =  f .  (V'  )  on 
ij  l  j  ~R  c 


Yr  *  Jf  -  y  =  vRt  . 


(79) 


where  n  is  the  unit  outward  normal  to  the  surface,  is 
the  relative  velocity  between  the  die  and  the  deforming  body 
and  t  is  the  unit  base  vector  along  the  die-workpiece  interface. 
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Let  v*  be  any  admissible  velocity  field  satisfying  the  incompres- 
sibility  condition  e*  =  0  and  the  velocity  boundary  conditions,  the 
virtual  work  principle  gives: 


a. .e?  .  dV  = 


dS 


■f 


dS  + 


(a.  .n.) V.  dS, 
ij  j  1 


(80) 


where  V*  =  v*  -  is  the  relative  velocity  between  the  deforming  mater¬ 
ial  and  the  die,  associated  with  the  admissible  velocity  v* . 

For  the  case  where  the  magnitude  of  the  frictional  stress,  f,  is 
specified  according  to  the  assumption  that  the  ratio  of  the  frictional 
stress  to  the  shear  yield  stress,  k,  remains  constant,  f  is  expressed  by 


f(V  )  =  -mk 


-R'  -  --  [VrT 


(81) 


where  m  is  the  friction  factor  (0  _<  m  1).  Eq.  (80),  together  with 
the  inequality 


mysfR  -  >  °- 


(82) 


according  to  Eq.  (81),  and  the  convexity  and  normality  conditions  of 
the  yield  surface. 


(a? . 
ij 


a.  .)e*  .  >  0, 
ij  ij  - 


(83) 


gives 


a*  £).  dv  -  [  F*v*  dS  ■ 

K  u  *j  Js  -  -  i s 

F  c 


f(V*)-y*dS>J  (aijnj)VidS, 


(84) 


+  t'(Vp)  *uu  dS 

S  ~R 
c 
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where  a*,  is  the  stress  associated  with  e*.  and  on  the  yield  surface, 
ij  ij 

Removing  the  constraint  of  incompress’bility  condition  of  an  admissible 
velocity  field  v*  (denoted  by  u)  by  introducing  the  Lagrange  multiplier 
A,  the  stationary  value  problem  for  the  finite-element  formulation  is 


now  stated  as 


TU  {j  51  dv  - 

j  v 


F’U  dS  - 


f(vR).yR  dv 


Aey  dVi  =  0. 


Eq.  (85)  is  the  one  used  in  [56].  Note  that  in  Eq.  (85)  the  relative 

velocity,  VR,  and  the  volumetric  strain  rate,  ev,  is  now  defined  by 
D 

YR  =  u  -  y  and  =  u^  ^ ,  respectively. 

In  order  to  have  the  stationary  value,  the  existence  and  continuity 
of  the  derivatives  of  4>  with  respect  to  u  is  required.  Therefore,  no 
converged  solution  can  be  expected  when  two  consecutive  nodal  point 
velocities  relative  to  the  die  become  zero  or  close  to  zero,  by  the 
direct  application  of  Eq.  (85).  In  order  to  overcome  this  difficulty, 
the  use  of  velocity-dependent  frictional  stress  is  proposed  as  the 
interface  condition.  It  is  expressed  by 


f  =  -mk'i  (z-)  tan’1  [■ 


l  l  D| 
^a;u  |- 


where  the  die  velocity,  u11,  is  equal  to  unity  and  a  is  a  constant 
which  is  several  orders  less  than  the  die  velocity  (say,  10'^)  so  that 
the  frictional  stress  given  by  Eq .  (86)  is  almost  identical  to  that 
of  Eq.  (81),  for  |VR|  slightly  larger  than  zero  (say,  |VR|  >  10'd) . 

With  the  interface  friction  given  by  Eq.  (86),  the  stationary 
condition  is  expressed  by,  instead  of  Eq.  (85), 
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3u 


_3_ 

3u 


ae  dV  + 


r 

Ac 

v 


dV 


F*u  dS 


(  f-dV  )  dS 
J0 


0. 


(87) 


For  the  proof,  note  that  the  following  inequality  is  used  instead  of 
Eq.  (82): 


/  f 

'•■'0 


|  V*  i 
:~R' 


f-dVr 


f 

Jo 


f(Vp) • (V* 


-  vR) 


(88) 


Discretization  formulations  of  the  last  term  in  Eq.  (87)  are  given  in 
the  following  subsection. 


B_.  Discretization  of  the  friction  term 

Denote  the  last  term  inside  the  brackets  in  Eq.  (87)  by 


[j0  f-ivR]  dS. 


(89) 


Let  Tj  and  r0  be  the  radial  positions  of  the  two  consecutive  nodal 

points  1  and  2,  whose  velocity  components  are  (u  ,u  )  and  (u  ,u  ), 

rl  Z1  r2  Z2 

respectively.  The  unit  base  vector,  t,  is  assigned  along  the  die- 

workpiece  interface  as  shown  in  Fig.  24.  The  relative  velocities,  VD 

~R1 

and  V  are  equal  to  the  nodal  point  velocities  in  the  r-direction, 

~*2 

respectively,  namely. 


~R 


1 


~R 


(90) 
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frt  f  Irt 


With  a  bilinear  velocity  distribution  for  a  quadrilateral  element, 
the  relative  velocity  along  the  interface  is  a  linear  function  of  r 
and  is  expressed  by 


For  the  frictional  stress  represented  by 


f  =  -(mk  1-1  tan'1  — \  t  . 


$f>  given  by  Eq.  (89)  can  be  written  as 


for  VR  £  0 ; 


for  VR  £  0  - 


(91) 


(92) 


(93) 


Substituting  Eq.  (91)  into  Eq.  (93)  and  performing  integration  and 
differentiation,  the  expressions  necessary  for  the  Newton-Raphson 
procedure  are  obtained  as  follows: 
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3.  COMPUTATIONAL  CONDITIONS 

The  ratio  of  the  outer  diameter,  the  inner  diameter,  and  the  thick¬ 
ness  of  a  ring  specimen  was  6:3:2.  The  computation  of  ring  compression 
was  performed  for  several  materials  under  various  values  of  m  in  Eq. 

(86)  for  interface  friction  (a  =  10‘5) .  The  materials  were  selected 
from  several  references  and  are  tabulated,  with  their  stress-strain 
relationships,  in  Table  7(a). 

Table  7(a) 


Materials  and  their  stress -strain  curves  used  for  computation 


No.  ;  Material  j 

Stress- 

1 

strain  relation  \ 

1 

Initial  yield 
stress 

Reference 

1  ! 

1  ;  Commercially  | 

5/V(1 

+ 

£  ,0.34  | 

0 . 3193' 

Yq  =  15,500  psi  : 

[48] 

i  pure  aluminum 

••> 

1 

(106.86  MN/m“) 

i 

2  O.F.H.C. 

5/V(1 

.  + 

£  .0.28  1 
0  3518 

Yq  =  42,000  psi 

[48] 

■  copper 

1 

(280.56  MN/iT) 

1 

3  |  A1 .1100F 

5/y0=O 

♦ 

0.768E) 

Yq  =  16,500  psi 

[56] 

1 

| 

!  (113.76  MN/m2) 

1 

j 

4  i  Annealed 

5/Vo-d 

+ 

£  ,0.3 

0 . 05205 

Yq  =  9, 100  psi 

[52] 

1  1100  Al. 

1 

1  (62.74  MN/iQ 

L_ 

The  stress-strain  curves  corresponding  to  the  materials  given  in  Table 
7(a)  are  shown  in  Fig.  25. 

The  friction  conditions  and  the  mesh  system  used  for  computation 
are  given  in  Table  7(b). 
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Table  7(b) 


Interface  friction  condition  and  mesh  system 


No.  Material 

Die-workpiece  interfacial 
friction  condition 

Mesh  system 

1  Commercially 
pure  aluminum 

m  =  0.2;  0.12  with  k  =  YQ//3 

m  =  0.12  with  k  =  a//3 

Fig.  7.3(b) 

2  O.F.H.C. 

copper 

m  =  0.6,  1.0  with  k  =  YQ//3 

m  =  1.0  with  k  =  a//3 

Fig.  7.3(b) 

3  A1.1100F 

m  =  0.5  with  k  =  YQ//3 

Fig.  7.3(a) 

4  Annealed  1100 

aluminum 

m  =  0.25;  1.0  with  k  =  a// 3 

Fig.  7.3(b) 

The  two  mesh  systems  used  in  calculation  are  shown  in  Figs.  26(a) 
and  (b) .  Both  systems  consist  of  96  elements  and  117  nodal  points  each. 
Fig.  26(a)  is  the  same  mesh  as  that  in  [56],  while  in  Fig.  26(b) 
smaller  elements  were  used  near  the  die-workpiece  interface  and  both 
corners  because  of  the  large  gradient  of  strain  rate  there  and  possible 
folding  formations  at  corners. 

The  non-steady-state  deformation  problems  in  ring  compression 
was  analyzed  in  a  step-by-step  manner  by  treating  it  quasi- linearly 
during  each  incremental  deformation.  The  reduction  in  height  at  each 
step  is  1%.  The  solution  of  velocity  fields  at  the  first  step  of  [56] 

(m  =  0.3)  was  used  as  an  initial  guess  for  all  the  friction  conditions. 
The  measure  of  the  convergence  is  represented  by  the  quantity  fl  Auff/llu/| 
where  Au  is  the  error  vector  found  during  iterations  and  the  Euclidean 

/n  2  2 

vector  norm  is  defined  by  ||u||  =  /  £  (u  +  u  ),  where  N  =  total  number 

1=1  i  i 
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Fig.  25  Stress-strain  relationships 
for  (1)  commercially  pure 
aluminum,  (2)  O.F.H.C.  copper, 
(3)  A1 .  1100-F,  and  (4)  Annealed 
al.  1100. 


Fig.  26  Mesh  systems  used  for  computation. 
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of  nodal  points.  The  solution  obtained  from  the  previous  step  was  then 
used  as  an  initial  guess  for  the  subsequent  step.  Normally,  8  to  30 
iterations  were  required  for  the  first  step  to  reach  the  accuracy  of 
Au  /  u  £  0.00005,  depending  on  the  difference  of  friction  factor,  m, 
between  the  present  case  and  initial  guess.  For  subsequent  steps, 
only  3  to  5  iterations  were  necessary  to  reach  the  same  accuracy. 

Some  additional  iterations  were  required  during  steps  where  inside  ov 
outside  folding  was  formed. 

4,  RESULTS  AND  DISCUSSION 

The  computed  results  are  shown  in  Fig.  27  in  terms  of  changes  of 
minimum  internal  diameter  as  functions  of  reduction  in  height.  The 
figure  revealed  that  the  geometrical  changes  differ  depending  upon  the 
selection  of  the  shear  yield  stress,  initial  or  current,  for  the  same 
m  value.  The  difference  becomes  significant  for  large  reductions  in 
height  for  both  m  =  0.12  and  m  =  1.0.  It  is  also  seen  that  the 
geometrical  changes  are  affected  by  the  material  property  as  indicated 
by  the  curves  of  m  =  1.0  for  copper  and  for  annealed  aluminum. 

A  comparison  of  the  present  numerical  solutions  with  those  by  the 
upper-bound  method  [49],  [50]  was  shown  in  Fig.  28.  For  low  friction, 
the  two  solutions  coincide  with  each  other  for  small  reduction  range, 
but  deviate  considerably  for  large  reductions.  It  is  interesting  to 
note  a  coincidence  of  the  two  solutions  for  m  =  0.12.  The  numerical 
solution  for  this  case  was  obtained  using  the  current  shear  yield  stress. 
For  high  friction  discrepancies  of  the  two  are  more  significant. 

Experimental  data  taken  from  the  literature  [47],  [50]  were  also 
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mum  internal  diameter, 


-  k  •  yA 


— -  Upper  bound  (Avitzur;  DePierre  et  al.) 


-  FEM  (k  ■=  3) 

-  FEM  (k  «  d//3) 

O,  •, A  Experiment 


Fig.  27  Changes  of  minimum  internal 
diameter  as  functions  of 
reduction  in  height 
(finite-element  solution). 


Fig.  28  Comparison  of  the  finite- 
element  solution  with  the 
upper-bound  solution  [49],  [50] 
Experimental  data:  o  -  copper 
[50],  a  -  copper,  •  -  pure 
aluminum  [48] . 
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plotted  in  the  figure.  Interpretation  of  this  particular  comparison 
for  low  friction  is  that  both  solutions  suggest  the  increase  of  friction 
as  reduction  in  height  increases.  However,  the  interpretation  of  experi¬ 
mental  results  for  high  friction  depends  on  the  theoretical  solutions. 
According  to  the  upper-bound  solution,  the  m-value  is  close  to  1.0  in 
the  beginning  of  compression,  and  the  friction  decreases  as  the  deforma¬ 
tion  increases.  On  the  other  hand,  the  present  numerical  solution  sug¬ 
gests  the  opposite,  namely,  the  m-value  is  about  0.6  in  the  beginning  and 
tends  to  increase  as  the  reduction  in  height  increases. 

In  Fig.  29  the  present  results  were  compared  with  another  upper- 
bound  solution  [51].  It  is  obvious  that  this  upper-bound  solution  is 
not  acceptable  even  for  practical  applications. 

Another  comparison  was  given  in  Fig.  30.  The  upper-bound  solution 
for  annealed  aluminum  by  Lee  and  Altan  [52]  is  in  good  agreement  with 
the  present  finite-element  solution.  Experimental  data  show  that  the 
m- values  for  this  experiment  remain  almost  constant,  and  m  =  1  for  the 
dry  condition  and  m  =  0.25  under  the  lubricated  condition.  The  load- 
displacement  relationships  computed  using  these  m-values  were  compared 
with  measurements  in  Fig.  31.  Although  the  theoretical  curves  pre¬ 
dicted  higher  load  values  than  the  experiment,  the  present  numerical 
solution  gives  closer  predictions  than  those  by  Lee  and  Altan.  Note 
that  the  displacement  readings  in  the  experiments  were  not  corrected 
for  elastic  deflection  of  the  testing  machine  as  reported  in  [52]. 

The  finite-element  method  not  only  provides  overall  quantities, 
such  as  geometrical  changes  of  specimens  and  load-displacement  curves, 
but  also  reveals  the  detail  of  the  deformation  characteristics.  Fig.  32 
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Fig.  29  Comparison  of  the  finite- 
element  solution  with  the 
upper-bound  solution  [51]. 
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Fig.  30  Comparison  of  the  finite- 
element  solution  with  the 
upper-bound  solution  [52]. 
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upper-bound  solution  from  [52] . 
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shows  grid  distortions  at  50%  reduction  in  height  under  various  inter¬ 
face  conditions.  For  low  friction,  the  outward  flow  dominates,  but 
with  increasing  friction,  the  neutral  flow  point  appears  along  the 
die-workpiece  interface,  resulting  in  both  inward  and  outward  barrel¬ 
ing.  Furthermore,  for  high  friction,  folding  occurs  at  the  inner  sur¬ 
face  as  well  as  at  the  outer  surface.  The  deformation  characteristics 
shown  in  Fig.  32  are  qualitatively  in  agreement  with  experimental 
observations  reported  in  the  literature  [47],  [53],  [54].  The  deforma¬ 
tion  patterns  during  compression  can  be  seen  in  a  series  of  velocity 
distributions  shown  in  Figs.  33-36.  The  observations  in  these  figures 
are:  (i)  For  low  friction  with  m  =  0.12,  both  inner  and  outer  surfaces 

form  the  outward  barrel  contours.  The  neutral  flow  surface  (the  surface 
along  which  the  radial  velocity  is  zero)  was  first  found  at  about  45% 
reduction  in  height  near  the  inner  corner  and  moved  outward  as  deforma¬ 
tion  continued;  (ii)  For  the  m-value  of  0.25,  the  neutral  flow  surface 
was  first  formed  at  the  inner  corner  as  early  as  at  3%  reduction  in 
height.  It  moved  outward  and  formed  a  complete  concave  surface  towards 
the  z-axis  at  20%  deformation.  This  surface  kept  moving  outwards,  but 
the  shape  changed  and  became  convex  at  50%  deformation.  The  motions 
of  this  neutral  surface  resulted  in  an  outward  barrel  contour  of  the 
inner  surface  before  20%  deformation,  and  a  flat  inner  surface  at  50% 
deformation;  (iii)  For  cases  of  high  friction  (m  =  0.6  and  1.0),  neutral 
flow  surfaces  moved  outward  before  20%  reduction  in  height,  and  then 
moved  in  the  opposite  direction  as  the  reduction  increases.  The  shapes 
of  these  neutral  surfaces  are  different:  for  m  =  0.6  a  convex  surface 
was  formed  before  35%  deformation,  and  gradually  changed  into  concave 
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R-COOROlNflTE 


Fig.  32  Grid  distortions  at  50%  reduction  in 

height  for  various  m- values,  (a)  m  -  0 
(pure  aluminum),  (b)  m  =  0.25  (annealed 
Al.  1100),  (c)  m  =  0.6  (copper), 

(d)  m  =  1.0  (copper). 


91 


COORDINATE 


Fig.  36  Velocity  distributions  at  various 
reductions  in  height  for  m  =  1.0 
(k  =  o//3) . 
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surface  as  deformation  continued,  while  for  m  =  1.0,  the  surface  re¬ 
mained  almost  flat  during  all  stages.  The  formation  of  folding  can  also 
be  observed  from  grid  patterns  in  both  cases;  comparison  shows,  at  50% 
deformation,  the  amount  of  folding  for  m  =  1.0  was  almost  twice  as  large 
as  that  for  m  =  0.6. 

The  effective  strain  distribution  corresponding  to  two  extreme 
cases  (m  =  0.12  and  1.0)  at  20%  and  50%  deformation  are  shown  in  Figs. 

37  and  38.  For  m  =  0.12,  the  strain  distribution  was  more  or  less  uni¬ 
form,  in  general,  except  for  some  degree  of  strain  concentration  at  out¬ 
side  corners;  the  minimum  value  was  found  near  the  inside  corner  of  the 
die-workpiece  interface.  For  m  =  1.0,  on  the  contrary,  it  revealed  large 
degrees  of  non-uniform  deformation.  The  strain  increases  gradually 
from  die-workpiece  interface  toward  the  center  axis  of  the  workpiece  and 
was  highly  concentrated  at  both  inside  and  outside  folding  surfaces. 

Figs.  39,  40,  and  41  show  the  die  pressure  distributions  at  various 
reductions  in  height  for  the  m-values  of  0.2,  0.6,  and  1.0,  respectively. 
The  neutral  flow  point,  neutral  flow  regions  (regions  where  the  relative 
velocity  is  equal  to  or  less  than  0.0005  for  the  unit  die  velocity) 
and  folding  zones  are  also  indicated  in  the  figures.  In  Fig.  39,  ini¬ 
tially  the  die  pressure  distribution  is  uniform  with  slight  increase 
toward  the  outer  periphery.  As  the  reduction  increases,  the  die  pres¬ 
sure  at  a  distance  about  one-third  of  the  contact  surface  from  the  inner 
periphery  increases  more  than  that  for  other  parts.  The  neutral  flow 
point  which  appeared  after  28%  reduction  moves  toward  this  point. 

A  so-called  "friction  hill"  was  then  formed  at  about  70%  reduction 
in  height,  in  which  the  peak  pressure  was  at  the  neutral  flow  point. 
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Fig.  37  Effective  strain  distributions  for  m  =  0.12 

(k  =  Yq//3)  at  20%  and  50%  reductions  in  height 


Fig.  38  Effective  strain  distributions  for  m  =  1.0 

(k  =  Yq/ /T)  at.  20%  and  50%  reductions  in  height 


Fig.  39  Pressure  distributions  along  the  die-workpiece 


interface  at  several  reductions  in  height  for 
m  =  0.2  (pure  aluminum,  k  =  Yq//3) . 


Sticking  region 
Folding  rone 


Pressure  distributions  along  the  die-workpiece 
interface  at  several  reductions  in  height 
for  m  =  0.6  (copper,  k  =  Y-//3) . 
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Fig.  41  Pressure  distribution  along  the  die-workpiece 
interface  at  several  reductions  in  height  for 
m  =  1.0  (copper,  k  =  6// 3). 


r 

With  higher  friction  (m  =  0.6),  the  initial  pressure  distribution  is 

I  uniform  (not  shown  in  the  figure),  but  then  the  pressure  becomes 

minimum  inside  the  neutral  flow  region,  as  seen  in  Fig.  40.  The  pres¬ 
sure  near  the  inner  edge  of  the  contact  surface  increases  gradually 
and  becomes  maximal  at  35%  deformation.  Then  the  neutral  flow  region 
began  to  deform  and  pressure  there  increases  at  the  same  time,  and 
finally,  a  "friction  hill"  was  formed  at  60%  deformation,  at  which  the 
neutral  flow  zone  became  a  neutral  point. 

For  still  higher  friction  (Fig.  41)  the  pressure  distributions 
and  their  changes  are  similar  to  those  for  m  =  0.6,  except  that  the 
neutral  flow  region  was  present  at  57%  reduction  in  height.  It  must 
be  mentioned  that  similar  distribution  patterns  were  obtained  for 
annealed  1100  aluminum  to  those  for  copper. 

5.  SUMMARY  AND  CONCLUSIONS 

The  formulation  of  the  rigid-plastic  finite-element  method  was 
modified  for  the  neutral  flow  point  problems  by  assuming  velocity 
dependent  frictional  stresses  along  the  die-workpiece  interface.  The 
method  was  applied  to  the  analysis  of  ring  compression  for  several 
materials  under  various  friction  conditions.  The  finite-element  solu¬ 
tions  were  compared  with  the  upper-bound  solutions  published  previously. 
It  was  shown  conclusively  that  the  solutions  by  Lee  and  Altan  [52], 
among  other  upper-bound  solutions,  are  closest  to  the  present  solutions 
in  terms  of  the  geometrical  changes  of  ring  specimens  and  the  load- 
displacement  relationships.  The  present  finite-element  solutions  re¬ 
vealed  the  detailed  mechanics  in  ring  compression,  particularly  with 
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respect  to  the  formation  of  neutral  flow  points,  neutral  flow  zones, 
and  free  surface  barreling.  It  was  found  that  the  dominant  factor  for 
deformation  characteristics  is  friction  at  the  interface.  The  materials 


property  has  a  significant  effect  on  the  geometrical  changes  but  not 
on  the  die  pressure  distribution  patterns. 

The  success  of  the  present  formulation  for  the  analysis  of  ring 
compression  promises  solutions  of  the  other  neutral  flow  point  problems, 
such  as  rolling  and  forging,  by  the  rigid-plastic  finite-element  method. 


Section  VI 


CLOSED-DIE  FORGING 

1.  INTRODUCTION 

Forging  in  closed  dies  is  an  important  operation  in  the  shaping 
of  metals  into  useful  objects.  The  process  involves  the  compression 
of  a  specimen  between  the  upper  and  lower  dies,  restriction  of  the 
metal  flows  to  fill  the  closed  cavity,  and  flash  formation  through  the 
gap  between  the  closing  dies.  In  order  to  completely  fill  the  die 
cavity,  more  material  than  is  apparently  needed  is  normally  used  in  the 
forging  process.  The  flash  due  to  this  excess  material  is  subsequently 
trimmed  from  the  forging. 

The  flash  formation  itself  can  also  create  resistance  to  the 
escape  of  materials.  It  will  build  up  a  very  high  value  of  the  forging 
load  and,  therefore,  force  the  material  to  fill  all  the  recesses  of  the 
die  cavity.  It  is  important,  then,  to  adjust  the  dimension  of  the  flash 
so  that  the  extrusion  of  the  materials  through  the  flash  opening  is 
more  restrictive  than  the  filling  of  the  die  cavity,  and  that  the  forg¬ 
ing  load  required  is  still  not  high  enough  to  break  the  dies.  Thus, 
the  aim  of  the  flash  design  is  to  provide  the  minimum  flash  required  to 
give  the  best  filling  of  the  die  cavity. 

The  prediction  of  the  load  in  closed-die  forging  is  usually  very 
difficult  to  obtain.  Therefore,  most  of  the  mathematical  solutions 
are  suitable  for  open-die  forging.  Prandtl  [57]  considered  an  ideal 
plastic  solid  in  plane-strain  between  two  rough  dies  and  obtained  solu¬ 
tions  for  the  stress  distribution  and  average  forging  pressure.  This 
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case  was  later  considered  by  Hill,  Lee,  and  Tupper  [58],  Green  [59], 
Alexander  [60],  and  Bishop  [61].  For  closed-die  forging  at  the  finish¬ 
ing  stage,  the  actual  plastic  deformations  were  found  mainly  at  the 
flash  portion,  with  exceptions  of  minor  adjustments  of  the  forging 
material  to  fill  the  die  cavity;  it  is  therefore  reasonable  to  apply 
the  results  of  the  open  die  forging  of  a  flat  plate  or  solid  disk  to 
the  studies  of  flash  deformation  in  the  closed-die  forging  process. 

The  approximate  solutions  for  forgings  of  flat  plate  or  solid  disk, 
based  on  a  simple  state  of  stress,  were  proposed  by  Seibel  [4],  Nadai 
[62],  Sachs  [63],  Schroeder  and  Webster  [64].  Solutions  of  closed-die 
forging,  based  on  the  slab  method,  were  given  by  Kobayashi  et  al.  [65] 
and  by  Kobayashi  and  Thomsen  [66],  in  which  they  assumed  that  all  plas¬ 
tic  deformation  are  restricted  to  the  flash  and  to  the  materials  of  the 
forging  from  which  the  flash  is  formed.  More  recent  applications  of 
the  slab  analysis  to  closed-die  forgings  were  made  by  Altan  and 
Fiorentino  [67]  and  Akgerman  and  Altan  [68].  The  upper-bound  methods 
based  on  certain  admissible  velocity  fields  were  also  applied  to  forging 
problems  by  Johnson  [69],  Kudo  [70],  and  several  authors  in  the  USSR 
[73] - [74] . 

In  addition  to  the  mathematical  analysis,  an  empirical  approach 
to  the  prediction  of  closed-die  forging  load  has  also  been  proposed 
by  Schey  [71],  in  which  the  formula  P  =  aA^C^.  is  used,  where  P  is  the 
forging  load,  a  is  the  average  flow  stress  of  the  metal,  At  is  the 
cross-section  area  of  the  forging  at  the  parting  line,  including  the 
flash,  and  is  a  constraint  factor  which  depends  on  the  complexity 
of  the  shapes  (see  Spies  [72]). 
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Experimental  studies  on  the  closed-die  forgings  can  also  be  found 
in  several  references  [75] -[78].  Gubkin  [75]  and  Semenov  [76]  con¬ 
ducted  experiments  to  determine  the  zones  of  various  kinds  of  plastic 
deformation  in  forgings.  They  investigated  the  effect  of  flash 
geometry  on  the  maximum  load  for  the  closed-die  forgings.  Kondratenko 
[77]  and  Mertens  [78]  conducted  experiments  for  the  studies  of  grid 
distortion  during  the  deformation  process  in  closed-die  forgings. 

Kasuza  et  al.  [79]  observed  the  material  flow  in  plane-strain  and  axi- 
symmetric  rim-web  type  forgings. 

The  present  studies  utilized  the  same  material  (pure  lead)  and 
dies  as  those  used  in  the  experiments  of  plane-strain  closed-die  forg¬ 
ings  by  Lyapunov  and  Kobayashi  [80] .  Their  experimental  observations 
were  compared  with  the  present  finite-element  solutions  which  include 
grid  patterns,  velocity  fields,  forging  loads,  and  flash  formations. 
Although  the  finite-element  analyses  of  open-die  forgings,  with  rela¬ 
tively  simple  geometries,  have  been  conducted  by  Zienkiewicz  and 
Godbole  [25]  and  Price  [22],  the  applications  of  these  techniques  to 
problems  of  closed-die  forgings  is  still  practically  nonexistent 
because  of  the  difficulties  in  dealing  with  the  complexities  of  die 
geometries,  the  large  rigid  portions,  and  the  rapid  changes  in  velocity 
fields.  The  present  study  revealed  the  possibility  of  the  applications 
of  rigid-plastic  finite-element  method  to  the  analysis  of  these 
extremely  complex  problems. 
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2.  EXPERIMENT  BY  LYAPUNOV  AND  KOBAYASHI  [80]:  DIE  AND  WORKPIECE 
GEOMETRIES 


A  schematic  diagram  of  the  apparatus  used  in  the  experiment  by 
Lyapunov  and  Kobayashi  [80]  is  shown  in  Fig.  42.  A  2-inch  (5.08  cm) 
long  lead  specimen  with  a  square  cross-section  measuring  1.3  in.  * 

1.3  in.  (3.30  cm  x  3.30  cm)  is  placed  inside  the  cavity  between  the 
upper  and  lower  dies.  The  cross-section  of  the  cavity,  when  two  dies 
are  closed,  consists  of  portions  of  a  tapered  rectangular  cross-section 
and  a  tapered  flange  extended  from  its  upper  face.  Specifications  of 
the  dimensions  of  the  upper  and  lower  dies  are  given  in  detail  in  Fig.  42. 

As  described  in  [80],  the  plane-strain  constraint  was  accomplished 
by  the  use  of  two  parallel  platens,  one  of  which  was  made  of  1-in. 
(2.54-cm)  thick  tublex  tempered  glass.  In  order  to  minimize  the  effect 
of  friction  between  the  specimens  and  platens,  each  end  surface  of  the 
specimens  was  lubricated  with  oil.  However,  the  contact  surfaces  of 
the  specimens  and  dies  were  cleaned  with  acetone  and  dried  before  each 
test . 

The  forging  operations  were  performed  by  a  160,000-lb  (72.7  T) 

Tinius  Olsen  testing  machine  with  a  slow  speed  of  0.05  in.  (1.27  mm) 
per  minute,  at  room  temperature.  A  series  of  photographs  of  grid 
patterns  were  taken  and  the  load-displacement  curves  were  recorded. 

The  grid  patterns  were  then  used  to  obtain  the  particle  velocities  by 
projecting  the  negatives  of  the  photographs  onto  a  sheet  of  tracing 
paper  of  10X  enlargement  and  connecting  the  grid  points  for  the  several 
incremental  stages  of  the  forgings. 
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Lower  die 


Dimens  ions : 


0.54  in.  (0.864  cm) 
0.375  in.  (0.953  cm) 
0.75  in.  (1.905  cm) 
0.782  in.  (1.986  cm) 
1.2  in.  (3.048  cm) 
1.207  in.  (3.066  cm) 

1.75  in.  (4.445  cm) 

1.76  in.  (4.471  cm) 


Upper  die 


Workpiece 

! 

I  I  1.3".(3.30cmJ 


Fig.  42  Schematic  diagram  of  plane-strain 
closed-die  forging. 


The  results  obtained  from  this  experiment,  including  grid  patterns, 
velocity  fields,  load-displacement  curves,  flash  dimensions,  and  the 
height  variations  of  the  specimen,  were  used  in  comparison  with  the 
present  results  of  rigid-plastic  finite-element  analysis. 

3,  ANALYSIS  TECHNIQUE  AND  COMPUTATION  CONDITIONS 

The  rigid-plastic  finite-element  formulation  with  modified  penalty 
constraint  (Eq.  (18))  was  used  in  the  present  analysis  of  closed-die 
forging  under  plane-strain  conditions.  The  problem  is  characterized 
by  a  large  near-rigid  zone  along  the  center  portion  of  the  cavity, 
especially  at  the  initial  die  filling  stage;  this  rigid  portion  is 
subjected  to  a  large  portion  of  force-free  prescribed  boundary  at  the 
upper  surface  of  the  workpiece  inside  the  die  cavity.  To  deal  with 
this  type  of  problem,  the  proposed  method  of  rigid  treatment  (Section 

II. 2)  was  used,  which  assigns  a  limited  value  of  the  effective  strain 

•  • 

rate  (say,  £q  =  10~'j  in  the  equation  of  stress-strain  rate  relation¬ 
ship  for  all  elements  lying  inside  the  rigid  zone,  i.e.,  for  elements 
with  the  computed  value  of  effective  strain  rate  £  less  than  the 
assigned  value  of  e^. 

The  specimen  used  in  the  present  analyses  is  made  of  pure  lead 
which  is  characterized  by  a  rigid-perfectly  plastic  (i.e.,  non-work¬ 
hardening)  material  behavior  with  the  constant  flow  stress  Yq  assumed 
to  be  2,500  psi  (17.236  MN/m  ).  Two  extreme  cases  of  friction  conditions 
were  considered:  the  first  assumed  perfect  lubrication  between  the 
interfaces  of  the  die  and  workpiece  so  that  no  frictional  force  exists; 
the  second  assumed  a  sticking  condition  which  states  that  once  the 
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material  touches  the  die,  the  material  is  completely  adhered  to  the 
die  and  any  relative  movement  between  the  material  particle  and  the 
die  is  therefore  prohibited. 

The  mesh  system  initially  consists  of  208  elements  interconnected 
at  238  nodal  points.  A  very  fine  mesh  was  used  for  elements  lying 
below  the  inner  corner  of  the  upper  die,  and  also  for  elements  along 
the  free  surface,  as  shown  in  Fig.  43.  The  non-steady-state  forging 
process  was  analyzed  by  a  step-by-step  manner  with  the  size  of  incre¬ 
ment  equal  to  1%  of  the  initial  height  of  the  specimen  (i.e.,  0.013  in. 
0.033  cm).  An  arbitrary  initial  guess  of  the  velocity  fields  was  pro¬ 
vided  for  the  first-step  iterations.  The  solution  obtained  from  the 
previous  step  was  then  used  as  an  initial  guess  for  the  next  step. 

As  deformation  continued,  the  material  along  the  free  surface  moved 
towards  the  side  surface  of  the  die,  and  subsequently  touched  the  die. 
For  steps  where  the  die-workpiece  touching  phenomenon  did  not  occur 
during  the  previous  step,  only  5  to  7  iterations  were  required  for  the 
solution  in  order  to  achieve  an  accuracy  of  II A u II / II u II  <_  0.00008.  It 
was  found,  however,  that  the  velocity  fields  changed  dramatically 
for  steps  when  the  touching  phenomenon  occurred,  and  the  number  of 
iterations  required  to  achieve  the  same  accuracy  increased  accordingly. 

After  the  side  cavity  of  the  upper  die  was  almost  completely 
filled,  the  material  inside  the  lower  die  began  to  move  sideward  to 
fill  the  lower  cavity.  The  flash  formation  started  when  the  side 
cavity  of  the  upper  die  was  completely  filled,  and  it  became  appreciabl 
once  the  material  touched  the  upper  corner  of  the  lower  die.  This 
happened  at  about  the  27th  step  of  the  deformation  in  the  present  cases 
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Fig.  43  Mesh  system  for  plane 
closed-die  forging. 


In  order  to  properly  deal  with  the  deformation  mechanics  of  the  flash 
formation,  three  more  horizontal  grid  lines  were  added  to  the  mesh 
system  at  this  stage.  The  iteration  process  was  than  continued  until 
the  central  cavity  was  almost  completely  filled. 

Usually,  it  is  very  difficult  to  specify  the  boundary  conditions 
of  the  nodal  point  which  lies  at  the  corners  of  the  dies.  In  the 
present  analyses,  however,  these  nodal  points  were  assumed  to  be  fixed 
at  the  corners.  The  assumption  is  quite  reasonable,  as  may  be  seen 
from  the  experimental  photographic  grid  patterns. 

The  non-steady-state  grid  distortions  were  made  by  continuously 
modifying  the  grid  patterns  according  to  the  obtained  velocity  distribu¬ 
tions  at  each  step.  As  a  result  of  this  adding  process,  the  elements 
may  happen  to  cross  each  other  at  some  stage,  especially  for  the  ele¬ 
ments  around  the  inner  corner  of  the  upper  die.  In  order  to  avoid 
these  cross  elements,  the  grid  pattern  was  carefully  checked  at  every 
step  and  the  cross  elements  were  then  readjusted.  Note  that  very  tiny 
elements  were  also  found  around  this  inner  corner  due  to  the  large 
relative  displacements  of  the  surrounding  nodal  points.  These  problems 
may  be  solved  by  a  more  elaborate  device  of  the  element  arrangement, 
or  by  using  a  mesh  system  with  mixed  type  elements,  but  this  aspect 
of  the  problem  remains  for  further  investigations. 

4.  RESULTS  AND  DISCUSSION 

The  computational  results  for  both  cases  of  perfect  lubrication 
and  sticking  friction  wetemade  in  terms  of  grid  distortion,  velocity 
distribution,  load-displacement,  flash  geometry,  and  axial  height 
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variations.  Experimental  observations  in  [80]  were  used  in  comparison 
with  the  present  results. 

Grid  distortion 

Fig.  44  shows  a  series  of  grid  patterns  obtained  from  the  computa¬ 
tional  results  at  the  6th,  14th,  26th,  30th,  and  32nd  steps,  respectively. 
Each  step  represents  a  0.013-in.  (0.033-cm)  displacement  of  the  upper 
die  (the  lower  die  was  fixed) .  Experimental  observations  of  the  grid 
distortions,  at  approximately  the  5th,  16th,  23rd,  and  33rd  steps, 
a*.e  also  shown  in  Fig.  44.  The  comparison  shows  that  the  computed 
results  are  very  similar  to  the  experimental  observations,  and  the 
experimental  frictional  condition  falls  in  between  these  two  extreme 
cases.  The  following  features  were  made  according  to  the  provided  grid 
patterns:  (i)  At  the  beginning  of  forging,  the  specimen  is  compressed 

between  the  upper  and  lower  dies  and  the  material  at  the  upper  die 
corners  flows  sideways  to  fill  the  side  cavity,  while  a  large  portion 
of  the  material  along  the  center  of  the  specimen  remains  almost 
stationary.  For  the  case  of  perfect  lubrication,  the  movements  of  the 
nodal  points  along  the  upper  die  surface  is  faster  than  those  nodal 
points  at  the  side  surface  because  of  no  frictional  constraints.  The 
corner  cavity  of  the  upper  die  was  completely  filled,  at  the  8th  step, 
before  the  nodal  points  touched  the  side  surface  of  the  die.  For  the 
cases  of  sticking  friction,  on  Hie  other  hand,  the  folding  phenomenon 
was  found  first  and  the  side  nodal  point  touched  the  die  subsequently. 

The  sticking  friction  condition  resulted  in  a  small  die  cavity  at  the 
upper  die  corner,  which  cannot  be  completely  filled,  as  shown  in 
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(a)  Grid  patterns.  Upper:  experiment 
(5th  step);  Lower:  FEM  (6th  step, 
dark  portion  indicates  rigid  elements) 


[»; 


Sticking 


Frictionless 


Fig.  44  (b)  Grid  patterns.  Upper:  experiment 

(16th  step);  Lower:  FOM  (14th  step, 
dark  portion  indicates  rigid  elements) 
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Fig.  44  (c)  Grid  patterns.  Upper:  experiment 

(23rd  step);  Lower:  FEM  (26th  step, 
dark  portion  indicates  rigid  elements). 
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Fig.  44(b).  (ii)  When  the  side  cavity  of  the  upper  die  is  essentially 
filled,  which  was  at  the  26th  to  27th  steps  for  both  cases,  the  material 
starts  to  flow  through  the  gap  between  the  upper  and  lower  dies,  and 
a  small  portion  of  the  flash  is  formed.  This  metal  flow  is  then 
restricted,  and  the  filling  process  of  the  lower  die  cavity  becomes 
more  significant.  Note  that  a  small  cavity  at  the  lower  die  corner 
was  still  found,  at  the  32nd  step,  in  both  cases.  (iii)  When  both 
upper  and  lower  side  cavities  are  almost  completely  filled,  which  was 
at  about  the  32nd  step,  the  process  becomes  similar  to  extrusion  in 
which  a  large  portion  of  the  materials  are  forced  to  flow  upward  to  fill 
the  central  cavity;  the  flash  portion  between  the  upper  and  lower  dies 
increases  at  the  same  time,  (iv)  At  the  final  stage,  the  central 
cavity  is  completely  filled  and  the  process  involves  only  flash  forma¬ 
tion.  The  analysis  of  this  stage  requires  a  very  fine  mesh  system 
inside  the  flash  portion;  and  this  has  not  been  done  in  the  present 
investigation. 

Rigid  regions 

The  rigid  elements,  i.e.,  elements  with  the  computed  value  of 

-1  - 

effective  strain  rate  less  than  10  (e  _<  0.1),  are  also  shown  in  Fig. 
44,  for  the  case  of  sticking  friction.  The  dark  portions  in  the  figure 
represent  these  rigid  elements.  For  earlier  stages  up  to  the  7th  step, 
a  large  portion  of  the  workpiece  remains  rigid,  and  the  deforming 
elements  were  found  mainly  at  the  side  portion  of  the  workpiece  under 
the  upper  die;  this  deforming  region  enlarges  and  moves  downward  as 
the  process  continues.  For  deformations  from  the  10th  to  26th  steps, 


the  original  deforming  elements  under  the  upper  die  surface  become 
rigid  gradually;  this  rigid  portion  enlarges  and  finally  occupies  the 
entire  upper  corner  at  the  26th  step,  as  shown  in  Fig.  44(c);  the 
deforming  region,  on  the  other  hand,  shifts  to  the  lower  portion  of  the 
workpiece  and  the  elements  at  the  lower  corner  start  to  deform,  as 
indicated  in  the  figures.  Then,  the  lower  cavity  is  completely  filled 
and  the  deforming  elements  at  the  lower  corner  become  rigid  again; 
the  deforming  region  is  now  mainly  at  the  central  portion  of  the 
workpiece,  as  shown  in  Fig.  44(e). 

Velocity  distribution 

Velocity  distributions  corresponding  to  grid  patterns  described 
above,  for  the  case  of  sticking  friction,  are  given  in  Fig.  45, 
together  with  the  experimental  measurement  of  the  flow  patterns  from 
[80].  Fig.  45  shows  that  the  computed  velocity  distributions  are  in 
good  agreement  with  the  experimental  observations.  Material  particles 
at  the  upper  corner  moves  sideways  in  the  earlier  stages  and  changes 
its  direction  downward  in  subsequent  stages  when  the  particle  velocity 
of  the  lower  corner  becomes  appreciable.  The  particles  at  the  central 
portion  under  the  inner  corner  of  the  upper  die  start  to  rotate  and 
move  upward  when  the  lower  die  cavity  is  almost  filled;  the  particle 
velocity  at  the  lower  corner  becomes  zero  again,  at  the  same  time. 
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The  forging  loads  obtained  from  computation  are  compared  with  the 
experimental  results  in  Fig.  46.  The  predicted  values  for  both  per¬ 
fect  lubrication  and  sticking  friction  cases  are  lower  than  the  experi¬ 
mental  values.  The  deviations  may  be  explained  by  the  excess  load  re¬ 
quired  in  the  experiment  in  order  to  overcome  the  frictional  force 
between  the  specimen  and  the  constraining  plates  at  both  ends;  it  may 
also  be  due  to  the  additional  load  required  in  squeezing  the  materials 
to  flow  into  the  small  cavities  between  the  constraint  plates  and  the 
die  and  in  maintaining  the  plane-strain  conditions  for  the  forging 
process . 


Axial  height  variations  and  the  flash  dimensions 
Fig.  47  gives  the  variations  of  the  axial  height  of  the  workpiece 
under  various  deformation  stages.  The  computational  results  are  in 
good  agreement  with  the  experimental  measurements.  Note  that  the 
height  remains  practically  unchanged  for  the  die  stroke  up  to  10  mm 
(approximately  the  31st  step);  after  this,  a  relatively  large  particle 
velocity  in  the  upward  direction  was  found  at  the  upper  half  of  the 


central  position,  which  resulted  in  an  increase  of  the  axial  height 
of  the  workpiece,  as  can  be  seen  from  the  figure.  In  Fig.  48,  the 


geometries  of  the  flash  are  shown  in  terms  of  gap  distance  between  the 
upper  aid  lower  dies  and  the  dimensions  of  the  flash  width.  The  compari¬ 
son  shows  that  the  computed  results  match  well  with  the  experimental 


measurements . 
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Fig.  47  Axial  height  of  the  specimen  as  a 
function  of  die  stroke. 
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Fig.  48  Comparison  of  the  flash  dimensions 


5.  CONCLUSION 


The  finite-element  analysis  of  closed-die  forging  is  an  extremely 
difficult  facet  in  the  study  of  metal-forming  processes  because  of  the 
difficulties  in  efficiently  dealing  with  the  complexity  of  the  die 
geometries,  the  large  rigid  portions,  and  the  rapid  changes  in  velocity 
field  during  deformations.  Among  these,  the  rigid  problems  can  be 
solved  without  difficulty  by  using  the  proposed  modification  of  the 
stress-strain  relationship  (Section  II. 2).  The  present  analysis  util¬ 
ized  a  very  fine  mesh  system  in  solving  the  geometric  complexities. 

Grid  patterns  were  examined  and  readjusted  at  every  step  to  avoid  the 
possible  cross  elements.  It  was  found  that,  when  the  nodal  points 
touched  the  die,  the  iteration  numbers  increased  dramatically  because 
of  the  rapid  changes  in  velocity  fields.  Although  further  improvements 
are  still  required,  the  present  investigation  reveals  the  possibility 
of  solving  these  complex  problems  by  the  rigid-plastic  finite-element 
method . 

The  detailed  deformation  mechanics  of  closed-die  forging  under 
plane-strain  conditions  were  obtained  and  the  computed  results  were 
compared  with  the  experimental  observations,  which  include  grid  patterns, 
forging  loads,  variations  in  axial  height,  and  flash  formations.  The 
comparison  shows  that  the  computed  results  are  in  good  agreement  with 
the  experimental  data,  and  therefore  justifies  the  reliability  of  the 
present  rigid-plastic  finite-element  method. 

The  continual  changing  of  grid  patterns  in  closed-die  forgings 
revealed  the  following  deformation  stages  in  detail:  compression 
filling,  and  extrusion.  The  plastic  deformation  portions  of  the 
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workpiece  at  each  stage  were  further  explained  in  terms  of  the  distribu¬ 
tions  of  rigid  elements  and  the  velocity  distributions.  The  detailed 
information  on  forging  loads  and  flash  geometries  were  also  furnished 
by  the  present  method.  This  information  is  extremely  important  in  the 
selection  of  a  suitable  geometry  for  the  workpiece  and  in  the  proper 
design  and  control  of  flash  dimensions.  It  is  of  interest  to  note 
that  although  the  present  analysis  was  made  under  the  assumption  of 
two  extreme  friction  conditions,  the  rigid-plastic  finite-element 
method  itself  can  also  incorporate  any  other  friction  representation. 
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Section  VII 


SUMMARY  AND  CONCLUSIONS 


The  rigid-plastic  finite-element  method  has  been  successfully 
applied  to  the  problems  of  axisymmetric  extrusion  and  drawing,  ring 
compression,  and  closed-die  forging.  The  formulation  of  the  minimum 
principle  and  methods  of  imbedding  the  incompressibility  constraint 
were  discussed  in  detail.  Particular  attention,  in  the  finite-element 
formulation,  was  given  to  the  penalty  constraint  term,  when  the  Gaussian 
quadrature  numerical  integration  technique  is  employed.  The  original 
penalty  function  formulation  requires  reduced  numerical  integration  in 
order  to  avoid  the  possible  over-constraint  of  the  incompressibility 
conditions.  This  over-constraint  problem  can  be  resolved,  however,  if 
the  modified  penalty  function  is  used.  The  present  investigation  also 
suggested  an  efficient  method  in  dealing  with  the  rigid  elements  in¬ 
volved  in  the  computation.  This  method  of  rigid  treatment,  combined 
with  the  modified  convergence  criteria,  provides  fast  convergence  of 
the  iterative  solutions,  thereby  reducing  the  computing  time. 

The  four  nodal  quadrilateral  elements  with  bilinear  velocity 
distributions  were  used  throughout  the  present  investigations.  The 
discretization  procedures  change  the  minimum  or  stationary  value 
problems  into  a  set  of  simultaneous  nonlinear  equations.  These  equa¬ 
tions  were  solved  by  Newton-Raphson’s  iterative  technique.  The  time 
required  in  computation  for  each  iteration  depends  on  the  machine  and 
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the  number  of  elements  and  their  arrangement.  The  number  of  itera¬ 
tions  required  to  obtain  the  desired  convergence  solution  depends  on 
the  problem  itself,  as  well  as  on  the  initial  guess  provided.  In  order 
to  reduce  the  computing  time  without  sacrificing  accuracy,  coarser 
elements  were  used  when  the  strain  rate  gradients  are  small,  whereas 
finer  elements  were  necessary  for  large  gradients.  In  general,  the 
required  computing  time  is  reasonable  and  the  accuracy  of  the  obtained 
solutions  is  satisfactory  in  comparison  with  the  experimental  results. 

Both  steady-state  (extrusion  and  drawing)  and  non-steady-state 
(ring  compression  and  forging)  deformation  problems  were  analyzed  by 
using  the  present  finite-element  techniques.  The  non-steady-state 
problems  were  solved  by  a  step-by-step  deformation  process  in  which 
the  grid  distortions  and  effective  strains  were  obtained  by  incre¬ 
mentally  adding  the  computed  velocities  and  effective  strain  rates  at 
every  step.  For  the  analysis  of  steady-state  problems,  instead  of 
using  the  step-by-step  analysis  technique,  the  flow  lines  and  corres¬ 
ponding  effective  strains  were  constructed  at  every  iteration,  the 
flow  stresses  changed  accordingly,  and  the  iterative  procedures  were 
repeated  until  the  desired  converging  solution  was  achieved.  Work¬ 
hardening  material  behavior,  as  well  as  the  interfacial  friction  condi¬ 
tions,  can  also  be  treated  efficiently. 

The  analysis  of  two-pass  operations  in  extrusion  and  drawing 
processes  was  made  by  the  present  finite-element  method,  which  reveals 
the  effects  of  the  previous  deformation  history  on  the  flow  patterns 
of  the  current  pass.  The  deformation  characteristics  during  the  second 
pass  were  generally  closer  to  those  of  non-work-hardening  materials. 
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in  comparison  with  first-pass  deformation,  for  both  extrusion  and 
drawing.  The  axial  stress  along  the  axis  of  extrusion  and  drawing  was 
found  to  be  greater  in  tension  for  the  second  pass  than  that  during 
the  first  pass,  suggesting  an  increasing  possibility  of  center-burst 
formation  in  multipass  processes. 

For  neutral  flow  problems  such  as  ring  compression  and  rolling, 
the  boundary  conditions  along  the  die-workpiece  interface  are  not  well 
prescribed.  The  difficulties  were  resolved  by  assuming  the  velocity 
dependent  interfacial  frictions.  The  corresponding  minimum  principle 
and  discretization  formulations  were  made  and  applied  to  the  problems 
of  ring  compression.  The  results  revealed  that  the  inner  diameter 
changes  of  the  ring  specimen  are  highly  dependent  on  the  friction  condi¬ 
tions.  Although  material  properties  appear  to  affect  the  geometry 
change,  it  is  found,  however,  that  their  effects  on  the  velocity  and 
die  pressure  patterns  are  insignificant.  A  series  of  velocity  patterns 
were  also  given,  indicating  the  detailed  material  flow  under  the  in¬ 
fluence  of  friction  conditions.  The  present  results  match  well  with 
the  upper-bound  solutions  for  small  frictions,  but  differ  greatly  for 
large  friction  conditions,  which  indicates  that  most  of  the  upper- 
bound  solutions  did  not  handle  properly  the  situation  of  severe  barreling 
and  the  formation  of  surface  folding  under  conditions  of  large  friction. 

The  closed-die  forging  problems  were  analyzed  by  further  modifica¬ 
tion  of  the  program  for  two  extreme  cases  of  perfect  lubrication  and 
sticking  friction.  The  series  of  grid  patterns  and  velocity  distribu¬ 
tions  were  presented  and  compared  with  the  experimental  observations. 
Results  indicate  that  the  forging  process  may  be  characterized  by  the 
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stages  of  compression,  filling,  and  extrusion.  Each  stage  has  its  own 
deformation  characteristics,  as  can  be  seen  from  the  figures  of  rigid 
portions  and  velocity  distributions.  The  forging  loads  and  flash 
geometry  obtained  by  the  present  analysis  are  very  useful  for  the 
selection  of  a  suitable  geometry  of  the  workpiece  and  for  proper 
design  and  control  of  the  forging  processes. 

The  successful  application  of  the  present  rigid-plastic  finite- 
element  method  reveals  a  new  avenue  in  the  course  of  analysis  of  vari¬ 
ous  metal-forming  processes.  Results  obtained  from  this  method  are 
extremely  valuable  in  preparing  the  workpiece,  designing  the  equipment, 
controlling  the  process,  and  preventing  product  defects.  Although  the 
present  work  was  conducted  for  rate-insensitive  materials  under  room 
temperature  operations,  the  modification  of  the  present  method  to 
incorporate  strain-rate  and  temperature  effects  is  possible.  This 
may  require  a  new  variational  formulation,  but  the  basic  principles 
and  the  finite-element  techniques  remain  the  same.  The  present  investi 
gation  employed  the  friction  law  of  constant  interfacial  friction 
factor,  but  it  is  also  possible  to  incorporate  any  other  friction  repre 
sentations.  Although  further  investigations  are  still  required  in 
order  to  solve  complex  problems  with  a  minimum  of  computing  time,  it 
is  appropriate  to  make  the  conclusion  that  the  present  rigid-plastic 
finite-element  method  is  indeed  an  efficient  and  accurate  method  and 
is,  therefore,  a  powerful  tool  for  industrial  applications. 
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